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Preface 


The Present thesis an out come of researches carried out by me in the field of 
“Self similarity solutions for spherical & cylindrical shock waves in magnetodynamics” 
under the supervision of Dr. Kishore K. Srivastava, Head of Mathematics Department, 
Bipin Bihari (P.G.) College, Jhansi, alEfiliated to Bundelkhand University, Jhansi is being 
submitted for the award of Ph. D. degree in Mathematics. The thesis has been divided into 
eight chapters, each chapter has been subdi vided into a number of articles. 

The first chapter is introductory. It ^ves in brief, an idea about shock waves, 
spherical and cylindrical shock wave. Similarity principle, equation of motion and jump 
condition behind the shocks in magnetogasdynamics, radiation phenomennon, concept of 
self gravitation, magnetohydrodynamics waves, MHD and flow stability, MHD and Sun, 
characteristic method and Witham’s rule. 

The second chapter is devoted to unify study of self similar model of exponential 
spherical, cylindrical & plane shock wave taking magnetic radiative heat flux into 
account where total energy of the wave is variable and atmosphere is uniform. 

The third chapter deals with self similar study of power driven isothermal flow 
behind cylindrical shock. A qualitative behaviour of gaseous mass has been investigated 
with the help of «}uation of motion taking monochromatic radiation with gravitational 
effect in non uniform atmosphere. A comparative study is made between flow variables. 
Numerical integration of flow variables is carried out through software matlab. 

The chapter four, an attempt is made for the irtudy of propagation of plane shock 
wave in magnetogasdynamics under isothermal condition where temperature gradient is 
zero and radiation effect have been taken out. Numerical integration of differential 
equation is carried out on DEC 1090, computer at IIT Kanpur, by R.K. G’s program. 

The chapter fifth deal with “Self similar solution of cylindrical shock wave in 
magnetogasdynamic”. Taking gravitati(Mial foiw und^ non uniform atmosphere with 
radiation heat flux into account. The radiation pressure and radiation energy have been 


Ignored we have considered the problem shock wave reduced by a point explosion in a 
gas. 

The chapter six an attempt is made for study of the “Analytical solution of 
spherical shock wave in a rotating gas with axial component of magnetic field”. Taking 
density as constant ahead the shock front. By using C.C.W. method analytical solution are 
obtained for shock velocity and shock strength for weak and strong magnetic field. For 
strong shock, also we have considered two cases i.e. when the magnetic field is strong 
and when weak i.e. non magnetic case. 

In seventh chapter. We discussed with self similar model of radiative shock wave 
m magnetohydrodynamics with magnetic effect. The model of self similar shock waves 
produce on account of an instantaneous release of energy in a non-uniform equilibrium 

conditions. The disturbance are headed by a shock of variable strength. 

In eighth chapter an attempt is made for similarity solution of cylindrical shock 
wave with radiation energy and material pressure in magnetohydrodynamics. A thorough 
analysis of self similar equation of dynamics taking material pressure & radiation flux 
into account have been studied. The radiation pressure and radiation energy have been 
also consider. The gas in the undisturbed field, is assume to be at rest. It is grey and 
opaque in non uniform atmosphere. 
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Chapter - 1 


INTRODUCTION 


The study of ‘Fluid Dynamics’ is directed to the behaviour of a fluid in motion. 
The liquid in motion, and gas states are referred to generally as fluids. Some of the 
notable examples of application are ; flow of water through pipes; motion of an aircrafts 
or a missiles in the atmosphere. The study also yields the methods and devices for the 
measurement of various parameters, e.g. the pressure and velocity in a fluid at rest or in 
motion. 


The first notable works on fluids appeared in the seventeenth century. In 1687, 
Newton in his book ‘Principal’ deals with the influx of a fluid, the resistance of a fluid 
and the resistance of projected bodies and the wave motion. Daniel Bernoulli, in 1738 
determined the relation between pressure and velocity and formalized this in his theorem. 
In 1743 his son John Bernoulli applied the momentum principle to infinitesimal elements. 
Euler in 1755, realized the significance of the Bemoullis’ work and used it as the basis, 
setting down the fundamental equations of motion and equation of continuity for an ideal 
inviscid fluid. An alternate form of these equations was given by Lagrange, in 1781 and 
1789. Equations including viscosity were derived by Navier in 1822 and Stokes in 1845. 

In 1858 Helmholtz published a paper on vortex motion and in 1868 another on 
free stream line potential flow. 

Near the end of the nineteenth century, fluid flow itself began to be extensively 
observed and investigated. The modMii trend in fluid dynamics is to investigate the flow 
of electrically conducting fluid at very high temperature. 

At a microscopic level, the three states of a given substance are different because 
of the difference in the intermolecular distance. In many cases at normal conditions the 
molecular distance of a fluid are less than the minutest of any physical dimensions of 
practical interest. As a result, we interred in the st^istical average properties and the 
behaviour of large numbers of molecules, and not in that of individual molecules. That is, 
macroscopic and not microscopic, fM’opearties are of interest. In fluid dynamics as 
individual molecules are not being considraod, the fluid can be r^arded as a continuous 
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medium and so the physical quantities such as mass, momentum etc. of the fluid 
contained in a very small volume are regarded as being spread uniformly throughout that 
volume. 

In dealing with the gases at very low pressure, as in the upper atmosphere, or at 
very high temperatures such as in a plasma, the continuum concept of fluid dynamics 
must be violated and the study that has to be based on the behaviour of individual 
molecules (i.e. on the macroscopic approach) [1]. 

In our studies we frequently refer to a ‘small element of fluid’ which is always 
supposed so large that it still contains a very large number of molecules, as fluid is a 
continuous medium. So when we take of an infinitesimal element of volume we always 
mean that which are physically infinite small. Such an element is called a fluid particle. 

In continuum dynamics we assume that the macroscopic fluid properties, for 
example mean density, mean pressure, vary continuously with (a) the size of element of 
fluid considered (b) the position in the fluid and (c) the time. In (a), the variation 
becomes negligible as the element is physically very small. Thus, fluid properties density, 
pressure and velocity are expressed as continuous functions of position and time only. On 
this basis, it is possible to establish equations governing the motion of a fluid, which are 
independent of the nature of the particle structure. So gases and liquids may be treated 
together. 

There are two distinct methods of specifying the flow field [2]. 

LAGRANGIAN METHOD: In this method, the flow variables (velocity, pressure and 
density) of a selected fluid element or particle are described. If r,, is the position of the 
center of mass of the fluid element at time to then the basis flow quantity in the 
Lagran^an description is the velocity V (r^ , t). This method is also referral is individual 
same face of change. 

EULERIAN METHOD: In this method, the flow quantities are described at all points of 
space occupied by the fluid at all times i.e. flow quantities are defined as the function of 
position in space (r) and time (t). The basic flow quantity is vector velocity V (r,t). 
This method corresponds to ‘local time rate of change’. 
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f (x,y,z,t+ 5t) -f (x,y,z,t) 

Hence the local time rate of change’ is given by. 


As the point P is fixed the local time differential operator - 

a 

by the moving fluid. 

Now let at time t + 5t the fluid particle, which was at 
originally, is in the position ( x + u 6t, y + v6t, z + w 5 t ) where 

velocity components at the position at time t. The corresponding char 

f{x+u5t,y+va,z+wa,t+a)~f{x,y,z,i) 

or f (r + 6r, t + 5t)-f (r, t) 
and rate of change is, 

f(r+5f, t + 5t)-f (r,t) 


IS not carried along 


This gives the ‘individual time rate of change’. As the point P is moving it gives the rate 
of change which is carried along by the moving fluid. 

Since, 

f = f(x,y,z,t) 
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If V [u, V, w] be the velocity of the fluid particle at P and 


dx dy dz 
— =u, — = v, — = w 
dt dt dt 

then, 

df _ af af af 

-___u+ v + -— w 
dt ax ay az 


f; 

dt 


+ 


4 + ^ 

dx dy dz 


{ui+vj + wk) 


or, 
df af 


H” 


(v . v) f 


dt at 

where v = ui + yj + wk 

The first term on the right hand side represents the local rate of change of f, and the 

second term the convective rate of change. 

If V is introduced for f in the above equation then the total derivative of velocity 
with respect to time is, 

_ dv av /_ _ 

In cartesion rectangular coordinates 

(v. V) = u ^ + V — + w — 

ax dy dz 

It is this term in the Euler relation for acceleration which is responsible for the 

non-linearities in the equations of motion of fluid dynamics. 




FUNDAMENTAL EQUATIONS 

In order to study the details of the fluid flow, we wish to find the density 
distribution, pressure distribution, velocity distribution, states of the fluid etc. at all point 
of space occupied by the fluid at m times. Hence a knowledge of three velocity 
components (u,v,w), the tenqjerature T, the pressure p, the density p, etc. of the fluid 
which are functions of position in space ( f) ie. (x,y,z) and the time t, is needed. We 
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obtain relations connecting these unknowns and which would explain the particular 

problem of the fluid motion. We call these relations as fundamental equations. 

The flow of a compressible, non-heat conducting and inviscid perfect fluid is 
described by the four physical variables pressure p, density p, vector velocity v, and 

temperature T. So there must be four fundamental equations to find these variables. These 
are [3], 

Equation of continuity 
+ div (pv) = 0 


Dv dv , . \ 

15i = a 

When the external force is negligible. 

Energy equation is 

DQ Dfi _ j_ ^ 

Dt Dt p Dt 

For the adiabatic flow 

Dh Dp ^ 

P — ^ = 0 

Dt Dt 

and equation of state for perfect gas, 

P = pRT 

when thCTe is no intermolecular attractions. 

Now we consider the behaviour of the fluid flow at high temperature, where some 
of the fluids become conducting and the interaction between motions of the conducting 
fluids and variations in electromagnetic fields may not be negligible. To describe this 
interaction phenomena of the flow field with the electromagnetic field we find the 
velocity vector v, pressure p, density p , temp^ature T, and the m^netic field vector 

H, which are the functions of the special coordinates (x,y,z) and the time t. To calculate 
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these unknowns we must find the relations which are the fundamental equations for a 
conducting compressible fluid. In these equations there are coupling terms between the 

electromagnetic and fluid dynamics phenomena. We shall derive and discuss these 
equations [4], [5], [6], 


EQUATION OF CONTINUITY 

This equation simply expresses the law of conservation of 
fluid entering a certain volume in space must be balanced by tha 
matter is neither created nor destroyed. 

We now formulate this principle mathematically. 

Let V be any arbitrary volume fixed in space, bounded 
containing a fluid of density p . The volume element 6V is sn 
regarded as constant through it. 


pdV . The rate of generation of 


The mass of the fluid with in the volume V is 


the fluid within the volume is. 


PoO=0, as volume is constant with respect to 


If the volume V otxupied by a nnyving fluid, the fluid enters V through parts of its 
boundary surface S and leaves through another part. 
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Let h be a unit outward normal vector drawn on the surface element dS. The 
normal velocity is n v . The total outward flux is 


Using Gauss’s theorem, 

= Iv ^ • (pv) dV 


The sum of the net outward convection of mass plus the rate of generation, 

of the fluid within the volume must be zero. 


Since this is true for arbitrary elementary volumes 


The Maxwell’s electromagnetic equations for a conducting medium are. 


where B 


Differentiating the fourth Maxwell’s equation 


Using the second Maxwell’s equation it becomes, 
div (curl H-45tj)=47t^ 
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Third Maxwell’s equation is 


If the charged particle moves with velocity v 
J = qv then. 


This is the equation of continuity for the electric charge moving under the effect 

of magnetic field which is similar to the equation (1 . 1). 


EQUATION OF STATE 

The electrodynamical state relation is a simple relation between the current 
density J , fields and fluid motion, or the ohm’s law is, 

J = a ^ + p V X h} + q V 

Where 0 is the electrical conductivity and qv is the curreht depending on the 
motion of the net charge q. 

The state of a compressible fluid is defined by the pressure p, the entropy S, the 
internal energy s, the absolute temperature T, the mass density p and the specific 
volume V. 

An equation of state is formed by expressing any one of the thermodynamic 
variables in terms of the otl^r two quantities mentioiMd above. 

For example, ' 
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For a perfect fluid (no intermolecular intractions) having negligible viscosity the 


equation is, 


where, R is the universal gas constant. 

Assuming the process to be adiabatic and isentropic we have. 


This is the calorie equation of the state of the medium, s is the internal energj 
mass. 

When a fluid element changes its state isentropically the state equation is, 


EQUATION OF MOTION 

This ^nation express the laws of conservation of momentum in the fluid. The 
obtain equations deriving this we proceed as following 
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We consider a fluid mass in motion at time t occupying a volume V and bounded 
by a surface S. Let §V be a small element of volume. If the fluid has the density p , the 
mass of the element is p6V and it moves with velocity v (r, t) 


The inter force on the element is p6V 


This force equals the rate of change 


of linear momentum of the element. For the whole fluid mass. 


The force on the body is the sum of body and pressure force. 

Consider a surface element 6S and that n be the unit outward norma! vector to 
5S . Let p is the fluid pressure at any point on the element 8 S. Then the force on it due to 
the fluid outside the particle is, 

P (- n) 6S (For pressure acts along inward) Hence the total pressure force is, 

-J|npdS 

using Gauss’s theorem, 

Fp =-/£ gradpdV 





■' ' . ‘ ; 



10 


Chapter-I 



The electric body force on it due to an electric field of intensit 

The current flowing through the element is i = J 5A. Where 
section of fluid element. 

Hence using Biot-Savart law, magnetic force in the element is. 


Hence the total body force Fb on the fluid mass is. 


The above equation is true for arbitrary elementary volumes and so. 


This is the general equation of motion for a conducting fluid. Where Eulerian derivative 
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So the modified form of equation of motion is. 


can be neglected as compare to curl H as, 


Then the magnetic force is, 


Since (using vector identity) 


So equations of motion can be written as. 


In summation convention form it can be written as. 
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where 


EQUATION FOR THE VARIATION OF THE MAGNETIC FIELD 

This equation expresses how the magnetic field vary or the time dependei 
the magnetic field. 

The Maxwell’s electromagnetic field equations for a conducting medium are, 
div E = 4710 


conduction current J and the displacement current 


By n^Iecting the tarm 


above equation reduces to. 
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Taking Curl of this, 


where, J = o E + 


since, 


We may write, 


where 


In summation convention form it can be writtoi as. 


This is the magnetic field equation for a conducting compressible fluid 
For a fluid at rest it reduces to the diflusion equation. 
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-n IS called magnetic diffUsivity. For an equation of this type the rate of decay of the 
magnetic field is very rapid. 

And when the conductivity is so large that the term ^ V' H can be neglected that 
behaviour ofthe magnetic field is given by, 


In this case the rate of decay of magnetic field is much less therefore this is the 
case of high conductivity. 


EQUATION OF ENERGY 

The law of conservation of energy leads to another equation of fluid flow the 
energy equation. The law of conservation of energy is equivalent to the first law of 
thermodynamics which states that if a small quantity of heat is added to a simple system it 
is used up in changing the internal energy ofthe system and in the external work done by 


Where dQ is the heat added per unit mass, de is the increase in internal 


per unit mass and pd 


is the amount of work done during the change of volume 


The heat dQ is obtained by heat conduction 


For a perfect gas, the molecules have very negligible volume and there are no 
mutual attractions between the individual molecule and hence no potential energy. 
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then energy per unit volume of fluid flow is the sum of the kinetic energy density 
due to trauslaton- motion of the molecules i p and the internal energy density pe . 
The equation of motion is reduced to 


Dt P+ qE + p J X H, neglecting the gravitational force. 

Taking the scalar product of the above equation and the velocity vector v we get 
for the kinetic energy. 


The internal energy of the perfect gas depends on its temperature T only 


The equation of continuity for a compressible fluid is. 


It is an experimental fact that heat flow across an element of fluid surface in 6t 
time is proportional to the gradient in the temperature i.e. 
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The heat added to the system per unit volume is then. 

p6Q = 6t V [kVT] 

Therefore as 6t 0 


So for the internal 


Hence the sum of kinetic and internal 


Now we find the ordinary energy equation of an electromagnetic field. When the 
fluid moves in an electric and magnetic field it experiences a force known as Lorentz 
force and experimentally its value is 


Where qE is the electric force per unit volume due to an electric field and 
(j X b) is the magnetic force par unit volume due to magn^c field. 
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The rate of work done by the Lorentz force F is 


^ a»v,^vxn; = u ( 1 . 04 ) 

Now we transform the right hand side of above equation (1.04) into a function of 
E,D,B and H with the help of Maxwell equations. 

Since, 


Using the vector identity, we have 
E Curl H = H. Curl E - div (l x h) 
with the aid of this result 


If the medium is isotropic then 


18 


Chapter-I 


SH 


0, 


dt 

We write, 
E. J = 

or 


dt 


~ div (e X h) - ~ 
471 4tc 


dt 


(h.b)+| (I.d) 


E. J = - — div (e X h) 

471 ^ ^ dt 


d (pH^+eE 


2\ 


8x 


Since a part of the power density is the rate of variation of (e E^ -{xH^)/87t:, this 
quantity must be the sum of energy densities in electric and magnetic fields. 


d feE^ 


dt 


8x 


+ 


pH 

Stc 


2\ 


= _ f-£- 

^ 471 : 


div (e h) + E. j 


Using (1.03) and (1.04) equation (1.02) becomes, 


JD 

Dt 




€ + 


div(vp)+ F.v + div(kVT) 


(1.05) 


= [- div (v p) + E. J+ div (k V x)] 
Adding equations (1.05) and (1.06) we get 


(1.06) 


J?. 

Dt 


.2^ 


E + 


+ 


dt 


eE^ pH 
8n 8x 


2A 


V . (v p)+ V. (kvx) - A V (e X h) 

47C ^ 


This is an equation for total energy per unit volume of the flow field when there is 
no heating by radiant energy. 

The terms on the left hand side of above equation give the rate of change of 
energy per unit volume in the field. The right hand side terms of the above equation have 
the form of divergence so that, on integrating through a certain volume, all of them 
reduce to surface integrals. 

The first term represents the work done by the pressure force, second term 
represents the heat loss through conduction fi-om the surface and the last term represents 
the out flow of electromagnetic energy through the surfece where E x H is the Poyting 
vector. 
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The above equation of energy can be written as. 


For adiabatic flow no heat is added, conducted or radiated from the flow field. 


The energy equation for the adiabatic flow field 


SHOCK WAVES AND THEIR EXISTENCE 

If a small disturbance (i.e. with an infinitesimal amplitude and small velocity) is 
created within a non-viscous isentropic compressible fluid it will propagate throughout 
the fluid as a wave motion and with the velocity of sound relative to the fluid without 
suffering any distortion. 

In this case we get the wave equations which shows that both the density and 
velocity variations follow the same wave patterns derived from the linearizing of the 
equations of fluid motion. The solution p (x,t) and u (x,t) both are the single valued 

function. The equations of fluid nrotion reduced to linear form as velocity is small so that the 
term ( v grad) v in Eula^s equation of motion may be ne^ected. 

If the assumptions of infinitesimal amplitudes and gradients are removed then the 
wave velocity will not be a constant and a simple wave velocity will distort as it 
propagates. 

In this case as velocity is not small we have to solve the complete non-linear 
equations of fluid motion. The method of characteristics can be applied to solve such 
hypatolic type equations. 

A simple analytic solution of these equations of fluid motion may be find such 
that the density is a function of velocity only. After some manipulation we get equations 
in u and />. General solutions of these ^piations are • ; 

u = fi [x-(u±a)t] 




p = f 2 [x- (u±a) t] 

Where ft and h are arbitrary functions. 

We consider the solution, 
u = ft [x-(u + a) t] 

This equation shows that the disturbance is propagated at an instantaneous 
velocity u +a, instantaneous because this velocity is a function of time. If the increment in 
velocity u is much smaller than speed of sound then the solution is the case of sound 
wave i.e. the linear case and the curve u = n (x-at) does not change its shape as the 

disturbance propagates. But if u is not small the shape of the wave is disturbed as is 
propagates. 

The leailtirg (fistcxlKHis in vdoaty distrBxiticm are shown in %ore. 


Fig. -6 Formation of Shock 


At a time ti the points of high velocity move faster to the right than those of lower 
velocity. Thus the creast has moved faster to the right than the trough and so the profile is 
disturb^. The compressive part of the wave whCTe the propagation velocity is a 
decreasing function of x, distorts to give a triple valued solution for u (x,t) which is 
physically impossible for longitudinal wave. 

Actually the difficulty is over come by the formation of a 'shock wave'. In such 
cases a large variation in pressure and density occurs in a very narrow region in which the 
flow variables change rapidly and the fluid no longer undergoes isentiopic changes. The 
thickness of this shock wave region is veiy small and so we may considered the shock 
wave as a surface of discontinuity for many practical problems of inviscid fluids [4], [7]. 
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JUMP CONDITIONS ACROSS A DISCONTINUITY 
A differential equation can be represented as 


Where H, F, K are functions of position and time. K contribute for sources or 
sinks. F is the flux per unit time and H is the distribution of some state of medium. 
Functions H, F and K are continuous and diffia-entiable. 

Let X, (t)and X2(t) are continuous and differentiable function of t and 
X, (t) < Xj (t) for every t. 


Integrating (1.7). 
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Now we know that Newton-Leibnitz formula is 


Using this formula (1.08) becomes 


If X] and X 2 are constants then 


the outgoing flux and F (x, ) 


Where F 


is the incoming flux. 


The above result (1.09) states tha t 


The time rate of change of total amount of stale of the medium in any section 


is equal to the difference of outgoing flux and 


incoming flux and some contribution of sources and sinks, 


If there is no source or sinks i.e. K = 0 and. 


then (1.09) represents the conservation of mass. 
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and equation (1.07) is in the conservation form 


The most general conservation form for differential equation is 


suppose mere is a discontinuity at x = X (t) and xi and X2 are chosen so that 
X, <X (t) < X2. Suppose U" and and their first derivatives are continuous in X (t) 
> X > xi and in X2 > X > X (t) and have finite limits as x -^►X(t) from below and above. 


Integrating (1.10), 


Also using the property of line integral. 


Now we use the Newton-Leibnitz formula on the left hand side 


Where U 


and right respectively. 


)+ 

d 

f n ^ 



1 n + 1 J 
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where 


On comparing (1.11) and (1.12) we get 


is bounded in each of the intervals separately the integral tends to zero in the 


limit as x, x" and x, ^ x^ . Therefore. 


This condition may also be written as 


where the brackets [ ] indicate the jump in quantity and S is the speed of propagation of 
the discontinuity or shock front. 


This give the location of the discontinuity. So if the conservation equation is 
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then the jump relation is [7] 


JUMP CONDITIONS IN ORDINARY FLUIDS 


Consider a stationary shock front (discontinuity) separating two regions i.e, 

- = 0 so the Jump condition (1.14) becomes. 


The equation representing the conservation of mass, momentum and energy across 

a surface of discontinuity are given below 


Equation of continuity is 


which is in the form (1 . 1 3). So the corresponding jump condition (1 . 1 5) is 


Equation of energy is 


(using «|uation of continuity) 
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which is in the conservation form, jump condition is 


using ( 1 . 16 ) a we have, 


( 1 . 16 )b 


Equation of energy is 


In conservation form it can be written as 


where h = s + 
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Jump condition is 


PaUji n; then the constant factor may be dropped in the above 


condition and we have 


The equations [ 1 . 1 6] are the Runkine Hugoniot jump conditions when the shock is 
at rest The subscript land 2 denote a variable ahead and behind of the shock front are the 
components of the unit normal to the shock and bracket denotes the difference of values 

in the two sides of the shock surface of the quantity enclosed. 


JUMP CONDirnON IN MAGNETOGASDYNAMICS 

In an electrically conducting fluids in the presence of magnetic field, discontinuity 
(i.e a shock wave) in flow variables can exists. The study of magnetohydrodynamic shock 
waves was begun in 1950 with the paper of F. de Hoffmann and Taylor [9] .since then 
continued interest inspired by astrophysics by flight at the outer edges of atmosphere etc. 
has produced many papers describing shock wave properties. The basic properties of 
magnetogasdynamic shock waves are determined by the conservation laws. The presence 
of magnetic field the relations connection the flow and field quantities on the two sides of 
the shock surface are as follow [10] , Equation of continuity 
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which is in 


the form (1.13) so the corresponding jump condition (1.15) is, n^ ] 


(1.17)a 


The equation of motion is. 


since 


Using these results equation of motion can be written as 


which is in the conservation form. So the jump condition is, 
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Also we know 


Jump condition is. 


(1.17)c 


Jump condition is 


(1.17)d 


The equation of energy is 


V p + p8 + 


Ohm’s law for a moving medium is 
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also we know 


Curl H 


(neglecting the displacement current) 


In the case of high conductivity it reduces to 


using this result energy equation can be written as 


p+ps + 


p + pe + 


Jump condition is 
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The equation (1.17) are the shock conditions in the frame of reference in which 
shock is at rest. The equations (1.17) are also called as the MHD Rankine-Hugoniot 
relation, which can be written as from equations (1 . 17) 


For H = 0 this relation gives the result for a non conducting fluid 


The usual situation is that the flow ahead of the shock is known and these 


conditions are used to determine the flow behind or to determine the flow quantities in 
terms of one of the flow quantity behind. 


RADIATION PHENOMENON 

At very high temperature radiation can be considered as a continuous emission of 
energy in the form of electromagnetic wave which propagate in vacuum with the velocity 
of light. This energy is called as radiant energy of thermal radiation. Whereas, according 
to "quantum theory" the radiant energy emitted or absorbed is not continuous permitting 
all possible values, as demanded by the wave theory but in a discrete quantified form as 
integral multiples of an elementary quantum of energy, photon or light quanta. The 
amount of energy in each quantum being given by the product hv where O is planck's 

constant and v and fr«iuency of the radiaticm. 

Thus the quantum theory proposes the particle characteristics of radiation while 
the classical theory the wave charactmstics, both being required to understand the 
complex behaviour of radiation. It has also been recognised that ractiation reveals itself in 
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ct. cne eiecmcal (radio) waves, infrared, visible, ultraviolet. X-rays 
and gamma rays. This theoiy of thermal radiation can be applied to understand the 
processes which take place in stellar media, to explain the observed luminosity of stars 
and nuclear explosions and also to high temperature fluid flow. 

Radiative transfer and radiative heat exchange have an influence on both the state 
and the motion of the fluid. This influence is caused by the fact that fluid loses or gains 
energy by emitting or absorbing heat. The state of the fluid can be described by the 
fundamental equations which, in the presence of radiation .t,. 


RADIANT ENERGY DENSITY Er 

The radiant energy density per unit mass of the fluid is given by, 


Where a is stefan-Boltzman constant 


RADIANT PRESSURE Pr 

According to the well known result of classical electrodynamics the pressure of a 
radiation field is equal to the one-third of the radiant enerey densitv i e 


The radiant energy density and radiant pressure become comparable with the 
energy density and pressure of the fluid only at extremely high temperature or extremely 
low gas densities. 
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RADIATION FLUX F, 

Th© net amount of radiant energy passing 
unit time called the radiant flux through the surface 


Where C is the velocity of light in vacuum. 

The amount of energy lost or, conversely, the energy released in the fluid as the 
interaction with the radiation is given by the divergence of the radiation flux. Thus a 
comparison of radiant eneigy flux and the flux of energr of the fluid can be used to characteize 
the importance of the radiant heat transfer in the medium. 

Let, q be the energy lost by radiation in a unit volume of the fluid per unit time. 

then. 


If the fluid radiates more energy then it absorbs, then the energy is lost in radiation 
and q>0; if more energy is absorbed than emitted, then the fluid is heated by radiation and 
the energy loss is negative, q<0. 

Neglecting the effect of viscosity and heat conductivity, the fundamental 
equations for the fluid-radiation systems can be writtai as follow. 

The continuity equation remains unchanged 


The equations of motion taking radiant pressure into account and neglecting the 
external force are 
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replacing q by the divergence of the flux Fr 


SPHERICAL AND CYLINDRICAL SHOCK WAVE 


Consider the propagation of a shock wave through a perfect gas of intensity 
(i.e. very strong) resulting from a strcMig exploaoo i.e. from the instantan«)us rdease of a large 
quantity of energy. 

When the energy is suddenly released in an infinitely concentrated form and 
distribution of density, pressure etc. depends only on the distance from some point then 
this is the case of spherical shock wave. 

When the energy is suddenly released along a line and distribution of all quantities 
is homogeneous in some direction and has complete axial symmetry about that direction, 
then this is the case of cylindrical shock. 

Since we consider the symmetrical flow (centrally or axially) so there are only two 
indqrendent variables, namely r and t. So the one dimensional fundamental equations 
governing the adiabatic flow are [5], [14] 
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Where j 0, 1, or 2 stands for plane, cylindrical or spherical case, respectively, and r is 
the radial distance from the centre in spherical case, the radial distance from the line of 
explosion, in the cylindrical case u is radial velocity and p , the density. 

Similanty pnnciple may be used to reduce these equations to ordinary differential 
equations. 


SIMILARITY PRINCIPLE AND SELF-SIMILAR GAS MOTION 

It IS not always possible to solve non-linear differential equations describing the 
physics of a motion or a process by using mathematical techniques. An approximate 

solution of such a problem can be obtained by solving a similar problem which is easier 
to solve. 

We consider the problem of one-dimensional adiabatic flows of a perfect gas with 
constant specific heat with either planer, cylindrical or spherical symmetry. The system of 
equations for flow of this type is given by (1.18) to (1.20). These gas dynamic equations 
contain five dimensional quantities p, p, u, r and t. The dimensions of three of which 

are independent; for example density and time, so the equations admit three independent 
similarity transformation groups of quantities. Dimensional anaivc;<= 
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me mnction p (r,t) can be written in the form p (r,t) 

7c(t) P (r/R), where 7t(t) is the scale of the pressure and R (t) is the length scale both 

depend on time in some manner and the dimensionless ratio P = P (r/R) is a function 

71 

of new dimensionless coordinate „ = r/R. Multiplying the variables P and Oy the 

scale function it(t) and R (t), we can obtain from the function P ( ti), independent of time 

the true pressure distribution. The other flow variables, density and velocity, are 
expressed similarly [15], 

It IS a natural question that what requirement must be satisfied by the conditions of 
a problem in order that the motion be self-similar. Dimensional analysis is used to answer 
this question. Since the dimensions of pressure and density contain the unit of mass at 
least one of the parameters in the problem must also contain a unit of mass. In many cases 
this is the constant initial density of the gas p„ which has the dimension MUK Let the 
parameter containing the unit of mass is a. It can be assumed that its dimensions are 
[a]=ML''r. The dimensions of the functions p, p, and u are, 
[p] - M LT , Ip] = ML ^ and [u]=LT“*, we can without any loss of generality, 

represent them in the form suggested by Sedov [16], 
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The problem of a strong explosion represents a typical example of a self-similar 

motion. This problem was formulated and solved by Sedov [16] and succeeded in finding 

an exact analytic solution to the equations of self-similar motion. The same problem was 

also considered by Stanyukovinch in his dissertation [17] and by Taylor [18] both of 

whom formulated the equations for the problem and obtained numerical and not analytic 
solutions. 

The parameters in the problem of a strong explosion are the initial density of the 
gas p, ~ ML^ and the energy of explosion e~ ML^ . The energy E is alwavs eoual tn 


gas ana as a result an energy integral appears in the 

problem. These two parameters can not be combined to yield scaler with dimensions of 
either length or time. Hence the motion will be self-similar, that is will be a function of a 
particular combination of the coordinate r (distance from the center of the explosion) and 
the time t. The initial pressure and sneed of sound n and n in +r.n -j? . 


L T . Hence the dimensionless quantity 


can serve as the similarity variable. The distribution of pressure, density, 
velocity can be expressed as functions of one dimensionless variable n . 
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In the case of spherical symmetry the effect of all masses on a parti 
mass) at a distance r from the centre of symmetry is equal to the force of attracti 
point mass placed at the centre of symmetry and having a mass m (r, t). 

The behaviour of gravitationally interacting gaseous masses forming 
given by the appropriate equations of motion ([ 1 2], [ 1 6]). 

The continuity equation for radial gas motion with spherical symmetiy is 


For the mviscid gas, momentum equatitm taking gravitational forces into account 
(according to Newton’s law of gravitation gravity g at r is g = can be written as 


Where G is the gravitaUonal constant and m (r,t) is the mas sag as with in the 
sphere of radius r. 

To determine mass (m) we use the equation. 


which expresses the fact that the mass of a spherical shell of radius r and thickness dr 
4^ p and the energy equation for adiabatic gas motion 


A very large no. of papers have been published by now in whi 
similar solution were obtained to explain the adiabatic unst^y flow 
gas and analysed for systems of partial differential equations encoi 
problems of astrophysics such as internal motion in stm. 
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EQUATIONS OF MOTION AND JUMP CONDITION IN MAGNETO 
HYDRODYNAMICS 

In an electrically conducting fluids in the presence of magnetic field discontinuity, 
i.e. Shock wave in flow variable can exit’s. The study of magnetohydrodynamic shock 
waves was begun in 1950 with the paper of F.de Hoffmann and Taylor . 

When electric current induced in the fluid, then their flow in the magnetic field 
produces mechanical forces which modify the motion magnetogas dynamic owes its 
peculiar interest and difficulty to this interaction between the field and the motion. Thus 
the equations of magnetogasdynamic are the ordinary and electromagnetic equations. We 
only take into account the interaction between the motion and the magnetic field, we have 
Ignored the Maxwell’s displacement currents. Through out the thesis we take p = 1 
because magnetic permeability p differs only slightly from unity which is unimportant. 

We also assume that the dissipative mechanism such as viscosity, thermal conductivity 
and electrical resistant are absent. 

Since the problems dealt in this thesis relate to magneto gas dynamic or magneto- 
radiative shocks. We refer in this article the relevant flow and fluid equation. 

p = pRT; 

(1.18) 


Dp du . pu 
~f;~+p~+j~=0; 
Dt or r 


Du 4.1. 4. ^ ^ t)h 

Dt p ^ p 5r ^ 


Dh 5h 5u hu ■ 

-_4.u-— 4-h — 4-0 = 0; 

Dt Sr dt t 


|('’P‘’)+u|:(Pp'’) =0: 


D 


Where ~ = ^ + u ~ 
Dt dt dr 


and u,p, ,h,r and t are the velocity, jK^essure, density magnetic field transverse to 
the flow, radial distance and time respectively, T is the ratio of specific heat at constant 
volume and constant pressure, J=0, 1 and 2 corresponding to plane. Cylindrical and 
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spherical respectively and u=0 for plane and u=l corresponds to cylindrical and 
spherical both. 

In the presence of magnetic field the relation connecting the flow the field 
quantities on the two sides of the shock surface are as follows (34), (35) where the 

velocity in front the shock wave is zero : 

h2(u-U2) = h,u 


^ A 2 (r-l)A p, 

Wiete snfe 1 and 2 conespond to the value of the quantities just ahead and just 
shock sur&ce and u is the dKdc velodly. 


Slow movement of a fluid may be characterized by streamlines indicative of 
regular, orderly motion, called laminar flow. As the flow velocity is increased or the 
space available made larger there is a point where any small disturbance is amplified and 
the flow breaks up into tuihulence. Then the velocity profile, and all the other flow 
characteristics change radically. This is especially familiar in pipe flow, in the transition 
from laminar Poiseuille flow to tinfeulent Fanning-Hagen flow. The criterion is the 
dimensionless Reynolds number R = VLp/ fi , which indicates the relative importance of 

inertial and viscous forces. For flow in rough pipes, the transition occurs around R = 

2500 

A magnetic field may be expected to stabilize a flow against the transition to 
tuihulent flow. The cases of the magnetic field in the direction of flow and normal to the 
flow have been studied, and the stabilizing effect observed. The critical R^moWs raimber 
IS increased as the magnetic Reynolds number is increased. Experiments have been made 
that show this general result, but they are difficult and quantitatively not very Conclusive 
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Another kind of instability occurs when there is a temperature gradient and a 
pressure gradient in the same direction such as occurs in the lower atmosphere or in a pot 
of liquid heated on a stove. If we imagine a small parcel of fluid from one level suddenly 
moved to a level where the pressure is different, it will expand or contract adiabatically 
(to a good approximation) until the pressures are equalized, and its temperature will 
change by a definite amount. The change in temperature as a function of difference of 
position defines the adiabatic lapse rate, to use meterorological terminology. If the actual 
lapse rate (temperature gradient) is greater, then if the parcel has risen, it will find itself 
hotter and less dense than its surroundings. Buoyancy will then encourage it to rise even 
higher. Should the parcel have sunk, then things are exactly reversed, and it will find 
itself cooler and denser, so it will be impelled to sink further. A parcel at 1 just arrived 
from 0 is hotter than its environment 1'. while a parcel at 2 is cooler than its environment 
at 2'. If, on the other hand, the lapse rate is less than adiabatic, then a displaced parcel will 
be encouraged to return to its original position. In the former case, we have convective 
instability, and in the latter, convectil'e stohilill". Indeed, the lower part of the 
atmosphere, in which the decrease in temperature with height is usually more rapid than 
the adiabatic decrease, and which is. therefore, convectively unstable, is called the 
troposphere ("sphere of turning"). The unusual opposite case is called an invasion. 
Convection is important not only in the troposphere and in pots, but also deep in the earth, 
and in the sun. where energy is struggling to get out [19-20], 

A magnetic field in the direction of the temperature and pressure gradient will 
hinder the transverse motion essential to convwjtion, and make the convective "cells" 
narrower and less efficient, reducing the rate of oiergy transfer. In the Sun. this is seen by 
the relative darkness of sunspots, where a vertical magnetic field of thousands of gauss 
reduces the efficiency of convection of heat from below, so the surface cools below the 
general level of the photosphere. 


MHDANDTHESUN 

It was noted above that the tmly conducting fluids available fir laboratory 
experiments are mercury and liquid sodium, both inconvenient for different reasons. It is 
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very difficult to reach large magnetic Reynolds numbers in laboratory experiments, and 
so to verify important theoretical results with any accuracy. There is another MHD 
laboratory available, however, and that is the Sun. ITie surface of the sun is a hot 
relatively dense plasma where the magnetic Reynolds number is very large, so the 
magnetic field is well and truly frozen into this fluid. However, we cannot change the 
experimental parameters, and we do not know what is going on below the level that we 
can see. so it is a less than perfect laboratory. Nevertheless, there are many interesting and 
varied phenomena that show the influence of MHD very well. 

The visible surface of the sun is the photosphere, a fuzzy surface that we can see a 
little ways into, as into a cloud. It is composed of 90% hydrogen. 10% helium, with all 
the other elements as traces. Here energy rising from below encounters a steep declining 
temperature gradient as the surface layers cool by radiation into space. Convective 
instability produces active convection, with cells ("granules") averaging 700 miles in 
diameter, hot plasma rising in the center and sinking at the periphery after cooling a little. 
We observe a temperature of about 6000K looking straight down, closer to 5000K at the 
limb where we look obliquely and not so deep. The plasma emits a black-body spectrum, 
acting as an almost perfect emitter. The plasma cools rapidly above this level, becoming 
less dense and more transparent. In its lower temperature of about 4500K. neutral atoms 
absorb their characteristic lines, creating the Fraunhofer lines in the solar spectrum. This 
region is callal the revming laya-, and extends some hundreds of miles vertically [21], 

In solar eclipses, the moon covers the photosphere and the lower part of the 
reversing layer, but the upper part of this region is seen in its bright red Ha radiation at 
656.3 nra. which dominates the flash spectrum observed at this instant, when all the solar 
Fraunhofer lines become bright emission lines. For this reason, this region is called the 
chromosphere, the lower part of the solar atmosphere. Above the reversing layer, the 
chromosphere becomes rarer and hotter, merging with the whitish corona above 13.000 
miles. This is a rare plasma, in which the particles have veiy high kinetic enagies. 
characteristic of lO^K, but there is no thermal equilibrium here, so, temperature is a 
doubtful concept. The magMtic Remolds mimba- is high throughout the solar 
atmosphere, and matter clings to the lines of forces sliding along them fi-eely.' 
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The Sun is almost perfectly spherical, subtending an angle of about 32* at mean 
distance. The equator is inclined to the ecliptic, and we see the north pole in September, 
the south pole in March, from our position 1.495977 x 101^ cm distant. At the top of the 
atmosphere, we receive 1.37 x 10^ erg/sW, which coiresponds to a total luminosity of 
3.8 X 103^ erg/so The average density of the sun is 1 .4 g/cm2, about the same as Jupiter's. 
Its total mass is 2 x 103^ g, which makes its surface gravity about 28 times stronger than 
earth's. It is cunous that its equatorial period of rotation is 25 days (27 days, as seen from 
earth that is revolving about the sun), but it rotates more and more slowly towards the 

poles. At a latitude of 45°, the rotational period is 28 days, and has been said to approach 
33 days at the poles. 

This difference in rotational periods is noted without much comment in astronomy 
texts, but It IS really an extraordinary thing. In an ordinary liquid body, viscous forces 
would transfer angular momentum to the polar resioms until th<>v a* 
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... =H.»uuna, regions, but are restricted to mid-latitudes. Each maximum 

begins at high latitudes, and approaches the equator as the cycle progresses. The bigger 
the sunspot, the longer it lives. Large sunspots can live for weeks, and appear to cross the 
sun's disc more than once, but most do not. An average sunspot is 1000 miles in diameter 
miles in diameter, while a really large one could hold the earth. 

There is a number of phenomena surrounding sunspots and associated with the 
chromosphere and corona that make magnetic fields manifest. There are spicules and 
faculae and flocculae. coronal streamers and rays, all apparently supported by magnetic 
fields. The most impressive and long-lived are the filaments or prominences, best seen in 
elevation at the limb of the sun in their reddish Ha light, great arches of luminous matter 
supported by the magnetic field. The field itself seems to be the longest-lived of all the 
spot phenomena. The form of the coronal streamers suggested a dipole field, but the sun 
has no strong general field, which cannot be stroneer than 


CHARACTEMSTIC BIETHOD 


Many physical problems lead to the formulation of a quasi-linear system of first 
orda- equations. Such «}uations are linear in first derivatives of the dependent variables, 
but the coefficient may be functions of the dependent variable when these equations 
describe wave motion, a good understanding of many of the issues can be developed fi^om 
the study of plane waves. Accordingly, we start with the case of two independent 
variables. The two variables often the time and space variable so we denote them by t and 


45 





Chapter-I 


system is 


where the matrices A, a and vector b may be functions of vj - v„ as well as x and t. 

In general, any one of the equations has different combination of / Stand 
du^. Idx for Uj. That is it couples information about the rate of change of the different Uj in 

different directions. One can not deduce information about the increments of all the uj for 
a step in any single direction but we are at liberty to manipulate the n equal to see whether 

this information can be obtained from some combination of them we therefore consider 

the linear combination. 


where the vector / is a function of x,t,u and investigate whether / can be chosen so that 
( 1 .28) takes the form 


If this is possible (1.29) provides a relation between the directional derivative of 
all the Uj in the single direction (or,/?), when this is the case, it will be valuable to 
introduce curves in the (x, t) plane defined by the vector field (a, ;0) . If x= x(n), t = T 

( 7 ) in the parametric representation, of a typical member of this family, the total 
derivation of uj on the curve is 


withmit loss of gmioality, we may take 
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and write (1.29) as 


The condition for (1.28) to be in 


and we may eliminate the mj to give 


This is a condition on the direction of the curve such a curve in said to be 
charactenstic and the corresponding equation (1.30) in said to be characteristic form. 


WITHAM'S RULE 


cnester nas studied the motion of a shock wave down a non uniform tube on the 
basis of a linearlized theory in which the changes in the tube area and the consequent 
changes in shock strength have been assumed to be small. In this linearlized theory, the 
solutions breaks down when the flow behind the shock wave is nearly sonic Chisnel and 
Witham [7] have shown that Chesters work could be simplified and extended one minor 
simplification is that where as the Chester worked with the full three dimensional 
equations and performed and averaging process in the course of his analysis. It is 
sufficient to work from the outset with the one dimensional formulation. Witham has 
shown that the motion of the shock can be found in a simple way without solving the 
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equations for the flow behind the shock in detail. This method can be given as the 
following rule. 

The relevant equations of motions are first written in characteristic fonn. Then the 
rule is to apply the differential relation which must be satisfied by the flow quantities 
along a characteristic to the flow quantities just behind the shock wave. Together with the 
shock relations which express these values in terms of the shock strength, the rule 
determines the change in the shock strength. In fact it is found to the accurate even when 

the total changes m the shock strength are not small although a full understanding of this 
fact is still lacking. 
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Chapter - II 

Propagation of Exponential Magnetoradiative 

Shock Waves 


INTRODUCTION 


In this present chapter the self similar model of exponential spherical, cylindrical 
and plane shock waves is studied, taking magnetic radiative heat flux into account where 
total energy of the wave is variable and atmosphere is uniform. Similarty method for 
shock waves have studied wang [1], Koch [2], Helliwell [3], Ray and Baneijee [4] and 
other have investigated the propagation of plane shock waves in optically thick and thin 
limit cases of gas. In detail Gusev [5] and Ranga Rao and Ramana [6] have studied the 
problem of unsteady self similar motion of a gas displaced by a piston according to an 
exponential law. Verma and Singh [7] and Singh & Srivastava [8] have considered the 
problems of spherical shock waves in an exponentially increasing medium under the law 
of uniform pressure. Srivastava [9] has studied the problem of magnetoradiative shock in 
a conducting plasma. 

In this chapter we discussed the strong shock waves in a exponential spherical, 
cylindrical and plane shock waves in a uniform atmosphere with magnetic radiative 
effects. The similarity solution have been developed when radiation heat is more 
important than the radiation pressure & radiation energy and opaque the shock to be 
transparent and isothermal. The total energy of waves as cube of shock radius. 

The shock simulation is time dependent, and we assume the density to increase 
according power law. A solar flare accompanied by enhanced coronal temperature 
possibly exceedii^ 4x16^ gives rise to an ej^icm of plasma at a velocity of 500-1500 
kms'^ in interplanetary space, and to a shock wave propagating outwards fi-om the sun and 
reaching orbit of the Earth. In this case the shock wave will have a r^on where density 
increases in the direction of propagation J.B. Singh and P.R. Vishwakarma [10] have 
obtain and extent exponential shock with radiative heat flux. V.K. Singh [II] has find 
solution in conducting medium, R.F. Chisnell [12] has otrtain and extent an analysis 
descriptions of concerning shock waves. Dhee:^ Bhardwaj [13] has obtain and formation 
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of shock waves in magneto gas dynamics flows. Self similar flow behind a shock wave in 
a gravitating or non gravitating gas with heat conduction and radiation heat flux, were 
obtained by J.P. Vishwarkarma & Arvind K.C. [14], The shock wave propagate in an 
uniform atmosphere which is assumed to be at rest. The similarity solution have been 
developed when the radiation pressure and radiation energy is present viscosity and the 
solar gravity have been neglected. 


EQUATIONS OF MOTION AND BOUNDARY CONDITIONS 
The equation of flow behind a spherical, cylindrical and plar 
j — ^^®sponding to plane, cylindrical & spherical and u = 0 fc 
cylindrical & spherical both. 


(2.06) 

Where u, p , P, and T are the velocity, density, pressure and temperature 
respectively T being the gas constant are y the ratio of specific heat while F denotes the 
heat flux. 

Assuming local thermodynamic equilibrium and taking Rosseland’s diffusion 
approximation, we have 
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Where — is the stefan Boltmann constant. C the velocity of light, p the mean free path 
of radiation is the function of density and temperature following wang [1] we take 


° ^ (2.08) 
Po, a, &p are being constants 

The surface of constant discontinuity moves with time according to an exponential 

law. 


r=Aexp(mt) m >0 ^2 09 ) 

and since we have the flow is self similar, the shock will also moves with time according 
to an exponential law 


where r is the radius of inner expanding surface, R 
are dimensional constant. The disturbance is headed by i 
The boundary conditions are 


(2.14) 

0 and I denote the regions immediately ahead and behind shock front 

v is the shock velocity M denotes mach number. 


SIMILARITY SOLUTIONS 

The similarity transformation for the problems under consideration are 
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" " ^ ( 2 . 20 ) 

the variable n assumes the value 1 at the shock and u on the inner expanding surface. 

This enables us to express the radius of inner expanding surface. 


or shock velocity 


By equation (2.07) 


Using (2.21) & (2.22) 
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with P = -2,a remaining arbitary (0<a < 2) and 

M - _ 4- • . 

2p(p+4) ^ ^^^'^siOEicss rEdiation parameter. 


The equation (2. 15) to (2.20) are used to reduce the equation (2.01) to (2.06) and 
(2.23) in new form by equation (2.23), we get 


since 


By equation (2. 17) 
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from equation (2.01) 


substituting these values 


p&u inEq. (a) 


Po G r|m+ uvpg G — + Gov — + 


Equation (2.02) 


ucfr 1 aP h oh wh 

-j'w , . 

^ dr p dr p & pr 

By equation (2. 16) 
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since 


By equation (2.20) 


Substituting these values equation (b) 


mv-qV + mvV+mv W' + mv 


mH'q+ 


qV +V+W’ + 


P' = (HH’q-uH^J 

Equation (2.03) 


By equation (2.20) 
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By equation (2.16) 


By equation (2.20) 


Substituting these values in equation (c) 


Vp 7 vm [-H''n+H]+ vV^ vH' 

V. — [-H ti+h]+ vv— V 
r r 

Eh'ti + ^ + VH +V Hp+uHV 


Po V H vV 


uHV = 0 


By equation (2.04) 
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By equation (2. 17) 


By equation (2.19) 


By equation (2.18) 


P + vP 


Po |mvP+vP 


61 


Chapter-n 



again by equation (2.18) 


Substituting these values in equation (d), we get 


PG T| VPG T} Q't] jQ 


The shock conditions (2. 1 1)-(2. 14) are transformed into the following forms (at p = ) 


By condition (2. 12) 
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By condition (2. 13) 
Pi =Pov' 

Po v"P(ti)= Po v" 
at T} = 1 

p (i)= 1 

Condition (2. 14) 


1 


1 


y'M' 


1 


Po V 


Po v' Q (n) 


at T] =1 






PoV 


r M' 


-1 


By condition (2. 15) 
hi=yM^h, 

7p7 V H (ti) = y . V 


M 


at T]= 1 


H (l) = y M 

Thus the transformed equations of motion are 


V' 




(2.31) 


(2.32) 


(2.33) 


NP 


j PV + (q-V)H^ - oV-dH'-G 


[g (q - v)' - P - qH^(l - V)(q - V)J 

o(oV+jv) 

q(q-V) 

P’ = [HH'q - u J - G |v (q - V)- v) 
h' = (i-v)|hv + ''^ 


(2.24) 

(2.25) 

(2.26) 
(2.27) 
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The appropriate transformed shock conditions are 


Where primes denotes differentiation with respect 'p 


NUMERICAL SOLUTION ANB RESULT 

For exhibiting the numerical solutions, it is convenient to write the flow variables 
in the non dimensional forms as 


The numerical integration was carried out through using software matlab. 
certain choice of parameter mid reproduced in gr^hical form and nature of field varia 
is illustrated through them. We have calculated our result for foUowina Ham 


Ul 


p 

p 

Pi 

" Pw ’ 

h 

H 
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j 


(Pressure Distribution) 


■' — r — " — ^"T"' 


(Density Distribution) 


Spherical 

Cylindrical 


1 2 3 4 5 6 7 8 9 10 11 


Spherical 

Cylindrical 


1 2 3 4 5 6 7 8 9 10 11 
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Spherical 

Cylindrical 


(Magnetic Distribution) 


spherical 

Cylindrical 


(Velocity Distribution) 
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Chapter -III 

Self Similar Power Driven Isothermal Flow 
Behind Cylinderical Shock in Monochromatic 
Radiation with Gravitational force 

INTRODUCTION 

A self similar theoretical model of power driven isothermal expansion in a non 
homogeneous medium m the presence of monochromatic radiation with gravitational 
force IS considered. The result discussed depends upon variation of the flow variables 

behind the shock which are displayed graphically. Gas is assume to be gray and opaque 
shock. 

The theory of self similar flow behind shock waves in a non homogeneous 
medium has been discussed by many authors Sakurai [1], cole [2], Kopal [3] etc. Singh 
[4] has discussed the problem of propagation of shock waves in a one dimensional plane 
model of a non-homogeneous medium. Also Verma & Vishwakarma [5] have considered 
the problem of propagation of plane Cylindrical and spherical magnetogasdynamics, 
shock waves in a decreasing density medium using Whithams [6] rule. Sedov [7] made it 
possible to analyze certain classes of self similar solutions in a number of problems with 
disturbed energy release with the help of method of theory of dimensionally Onkar Nath 
[8] presented a theoretical model of cylindrical magnetogasdynamifccs shock waves 
under the action of monochromatic radiation in non-uniform steller atmosphere following 
work of Khudryautsev [9] on self similar problem of a motion of a gas under the action of 
monochromatic radiation. Nath [10], Manoj Pandey and V.D. Sharma [1 1] have discussed 
similanty analysis and find exist solution magnett^asdynamic equation further studied a 
model of cylindrical shock waves in a non-uniform rotating atmosphere under the action 

of monochromatic radiation, where he omitted magnetic field effect. 

We have studied the propagation of cylindrical shock wave in a 
magnetogasdynamics rotating non uniform atmosptere in the presence of monochromatic 
radiation and gravitation. The shock is assumed to be propagating in a conducting 
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medium at rest with density varing as {~2<p<Q). The magnetic field distribution 
varies as r“ (a <0) and is directed tangential to the advancing shock front. The radiation 
flux moves through the gas with a constant intensity in the direction opposite to that of 
the propagation of shock wave. Further, the rotating gas does not radiate itself and energy 

in absorbed only behind the shock wave. The radiation pressure and energy are very less 
hence neglected. 

For isothermal expansion of the plasma the internal heat generation per unit mass 
is identical for all elements of the fluid behind surface of discontinuity and is only a 
function of time so that for isothermal case. 


The medium in which expansion takes place is heterogeneous. The shock radius is 
a function of time. Thus density and shock radius have been assumed to obey the 
following power law. 


Where pg is the density of the ambient medium and R is the shock radius. Here A, 
B, a and P are constants. 


EQUATIONS OF MOTION AND BOUNDARY CONDITIONS 

In accordance with the above assumption the motion of an inviscid perfect gas in a 
magnetogasdynamics rotating non-uniform medium in prince of monochromatic 
radiation and gravitation can be described by the following system of differential equation 
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Treating the cylindrical shock front heading the expanding plasma as a surface of 
discontinuity. The jump conditions may written Whitam [6] as 
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Where suffixes 1 denotes state jump behind shock surface and ‘0’ der 

variables jump ahead of the shock surface and the dot’s stands for di 
respect to time. 

The Alfi^en mach number and usual mach number are defined as 


respectively where R 


is the speed of shock. 


The absorption co-efficient k, is considered as 
k = k, p" P” j^ r* t‘ 

where the dimension of constant ko is given by 


moreover the dimensions less constants Jq, Pq , Po are related as 


under the equilibrium condition, we 


SIMILARITY SOLUTIONS 

Let us consider the solution of the equation in the form 
u=RU(ti) 


where q is a non dimensional parameter d^ned to be 
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and thus p =1 at shock front r = R 


SOLUTIONS OF EQUATIONS OF MOTION 

Equations (3.01) to (3.07) may also be transformed with the help 
dimensional variables given in (3.22) - (3.30) 


By equation (3.23) 


By equation (ii) 
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By equation (3.23) 


By condition (3.22) 
u = R U (ti) 


substituting these values in equation (3.01) 


Po gRU' 


By equation (3.02) 


By equation (iii) 
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By equation (3.25) 


By equation (3.22) 
u = R U (ti) 


substituting these values in equation (3.02) we get 


H't} +RU 


■Po RH RU 


RH RIT 


H + H'(U-ti)+HU' + 


By equation (3.03) 


By equation (3.22) 
u = R U (ti) 
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By equation (3.24) 
P = p, X (rj) 


By equation (3.25) 


By equation (3.26) 
V = RV irj) 


substituting these values in equation (3.03), we get 
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Urj + U 


Po & nR 


By equation (3.04), we have 


By equation (3.26) 
V = R V (ti) 


again equation (3.26) 
v = RV(q) 
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substituting these values in equation (a) 


RV rjR^-tjR 


+RUt]R 


+ RV.RU = 0 


+ -nUV' + UV = 0 


V + q V 


By equation (3.05) 


Substituting putting these values in equation (3.05) 


2n g R^.tj R 
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W' = 


' (2 + p)g (2 + p)w 

n rj 


By equation (3.06) 


de de P 
— + u 
dt 


dr [St 


dp dp 
dr 


1 a V 

+ — ^ O’") = 0 

pr Sr 


de de 
+ u 


dp dp 
— + M — 
dr 


dt dr p^ ]_ dt 
By equation (3.29) 
e = R' E (q) 

— = 2RRE + R^ E'-^ 

5t a 

= 2Rfe^ :^E+ R" E' 
P R 


p dr pr 


v'Ry 


'2(13-1) 


a ~ R 

and 

dr dr 

dt E'R' 


E-qE' 


a 

Pr 


R 


ap Sp 

+u — 
dt dr 


^ X 
g' R 


-g'Tl + g 


a 

P 


j = Jo j(q) 

= '^0'^ 
dr R 

Substituting these values in equation (b) 


(p-i) 


E + qE' 


+ UE 


|g + (U-q)g' 


+ 




[J+XJ'] 


By equation (3.07) 


(3.35) 


(b) 


(3.36) 
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n+m-l 


where a, = k 


The boundary conditions are by equation (3.09) & (3.14) 


By condition (3.20) & (3.23) 


By condition (3.1 1) & (3.25) 


79 



Chapter-ni 


where M 


By condition (3.12) & (3.22) 


By condition (3.13) & (3.26) 


(r+i) 

By condition (3.14) & (3.29) 
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By condition (3.15) 


By condition (3.28) 
= Jo J (n) 
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RESULT AND DISCUSSION 

The set of differential equation (3.31) - (3.37) have been integrated numerically 
with the help of boundary conditions (3.38) - (3.45) by well known Runga-kutta method 

for a, = 2 and 2.5 p = 1.5 & -2, y=4/3, 7/5 ; = 5, = 10, a = „ = _y 2 

(y + l) 

M=3/2 ^ = 0, 5 = we have plotted the graphs showing the variations of various flow 
parameters with distance for different values of y, p and a, in pesence and absence of 
gravitational field and in presence and absence of radiation. This helped us to study the 
importance of gravitation and radiation respectively on flow parameters. 

From Fig. (3.11) to Fig. (3.18) we observe that the radial component of velocity, 
magnetic field, energy, rotational velocity, radiation flux decrease as we go towards the 
center of the explosion, while density and pressure increase as we go towards the center 
of the shock. In Fig. (3.18), it is surprising to note that the mass decreases for the case of 
y=7/5 for both presence and absence of gravitational field, it remains uniform for the 
case of y = 4/3 . In absence of gravitational field decrease in the radial velocity, magnetic 
field and energy is more prominent. In the presence of gravitational field we see that 
decrease in rotational velocity, radiation flux is more prominent while the increase in the 
value of density and pressure is more prominent. , , ■ " ■ 

To draw a comparison between gravitational effects vis-a-vis rotational effect, we 
have also observed the variations of flow parameters in the absence of rotation in fig. 
(3.19) to fig. (3.25) and compared it with the earlier drawn graphs for the absence of 
gravitation. We find that radial velocity, magnetic field, energy decrease more in presence 
of rotational velocity while radiation flux decreases more rapidly in the absence of 
rotational velocity . As expected, in the absence more rapidly in the absence of rotational 
velocity. As expected, in the absence of rotation, the value of pressure and density 
increase more as we go towards the center of the shock. However, the peculiar 
phenomenon to be observed is that the variation of mass remains uniform in case of 
7 = 4/3, a, =2 and 3 = 2 irrespective of the presence or absence of the rotational 
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effect. We further note that the gravitational effect is important for the propagation of 
shock waves in the present problem. 


We have calculated our result in suitable form 




0-8 0.85 0.9 0.95 1 

— p. 

fj 


Variation of Radial velocity with Distance series 1&2 shows 
with Gravitation (r=7/5. 4«) & series 3&4 without 
Gravitation (y=7/5, 4/3) 

Fig. 3.01 
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Variation of Magnetic fietct with Distance series 1&2 shows with 
Gravitation (7=7/5. 4/3) & series 3&4 without Gravitation (7=7/5, m ) 
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Variation of Magnetic field with Distance series 1 &2 shows with 
Gravi^ion (7=7/5, 4 / 3 ) & series 3&4 without Gravitation (7=7/5, 4 / 3 ) 
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Variation of Pressure with Distance series 1 &2 shows with 
Gravitation (y=7/5, 4/3) & series 3&4 without 
Gravitation ( 7 = 7 / 5 , 4 / 3 ) 
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Variation of Radiation flux with Distance series 1 &2 shows with 
Gravitation (y=7/5, 4/3) & series 3&4 without 
Gravitation (y=7/5, 4/3) 
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Variation of Mass with Distance series 1 &2 shows with Gravi 
tation (y=7/5. 4/3) & series 3&4 without 
Gravitation {y=7/5, 4/3) 
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Chapter - IV 

Pfopagation of plane shock wave in 
magnetogasdynamics 


INTRODUCTION 


iiie peeeoraena associated with heat transfer 
complex whether the fluid is at rest or in motion, study 
the piston problem with thermal radiation for 
similar method of Sedov [2], Helliwel [3] took 
with radiation heat transfer for general optically and 
the self similar solution for spherical blast 
approximation under the assumption that there 
symmetry laumnch and Probrteim [5] have studied the radiation 
an spherical medium several other author 


in rotating fluids are extremely 
or unsteady. Wang [1] considered 
one dimensional unstudy shocks using the 
a more general case of piston problem 
transparent limit. Elliot [4] discussed 
waves in air using the Rossland diffusion 
is no effect of heat flux at the centre of 
effects of shock waves in 
e.g. Summer [6]. Magnetic field also be taken 
into account. The unsteady model of Roche Consits of a gas distributed with spherical 
symmetry around a nucleous of large mass. It is assumed that the gravitating effect of the 
gas n^lected compared with the attraction of the heavy nucleus Rosenau [7] [8]. 

In this chapter we consider propagation of plane shock wave under isothermal 
condition. In the isothermal condition temperature gradient becomes zero behind the 
shock and radiation effect are already implicity involved. An idealized magnetic field is 
considered for only a portion of sphere including the origin i.e. the noint of exnlnrion 
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EQUATION OF MOTION FOR ADIABATIC FLOW 


where u, p, p, h and fare the velocity, pressure, density, magnetic field and radiative heat 

flux at a radial distance r from the center of the core at time t, m denotes the constant 

mass of the core, G be the gravitational constant, the magnetic permeability of the 
medium has taken to unity. 

The equation of state for ideal gas is given by 

^ = (4.05) 

Where r is the gas constant. 

Also assuming local thermodynamics equilibrium and taken Rossland’s diffusion 
approximation, we have 


Where (1/4) ac is the Stefan - Baltz man constant, C is the velocity of light and p the 

mean free path of radiation, is a function of density and absolute temperature T, following 
WangK.C. [1] 


where , a and P being constant, the total energy of the wave is dependent on time as 


Where Eo being the constant in the front of shock wave, we assume that 
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where A is constant and R denotes the radius of the shock surface ahead the shock, the 
magnetic field distribution is assumed to be as 


At the equilibrium state the other flow variable ahead the shock are 


where C, W and K are constant the disturbance is heated by an isotheraial shock and 
Ranking Hugoniot jump conditions are. 


2 “ ( 4 . 16 ) 

Where suffix 2 denote the flow variable just behind the and V be the velocity of the shock 
and is given by 


SOLUTION OF THE EQUATION OF MOTION 


In order to reduce the equation of flow to ordinaiy differential equation, we now 
introduce the following similarity transformation Sedov [2] ^ ^ ^ ^ ^ ^ ^ 
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We see the solution of equation (4.04 - 4. 19) in the form 


Where p is constant that assume the value 1 at the shock front hence we have 


By use of equation (4. 1 8) in equation (4.07) with the add of equation (4.05) we obtain 


Equation (4.01-4.04) and equation 6 are then transformed with the relation (4.18) and 
(4.20) to the form. 

du [Q(v- 5)] , 2{v5+l}->-^l}^-5)+(v+5).; 
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iP = p 1 _ Q 1 

dn II I v-5 LD“-p-"p^r:ij (;rr^ ^ 


1 dv 


m 

dri 


N 


(v-5) 


v(3-o) dv 
2tj dij 


dQ YP 


dri (r-l)ii 


(v-5)(2+Q) V ■ 

^ {®(t 1)+3y-1}+2 


Q 


^ (2-«>) + *(r-5) 


where L= 


4cap, 


(agm)“‘ 


3 j 13+4 

is a dimensionless radiation parameter. 


(4.24) 


(4.25) 


(4.26) 


(4.27) 


EQUATION OF MOTION FOR ISOTHERMAL FLOW 

The equation of continuity, momentum and magnetic field are the same as in the 

adiabatic flow but the energy equation becomes 
T = Constant 


dr 


0 


(4.28) 


By use of equation 4.05, we have 


P2 P2 


(4.29) 


By substitution of relation (4.18), (4.20) and (4.22) in «juation (4.01-4 03) and (4.27) the 
following set of equation are obtained. 

.jv __ 1 p0- v- 5)+ (v-5Xl + (o)-v]+N^X[v(3-a))-(l + G))(v-5)1 

~ {q(v-6)^D-N^-p| : 


4- 


(v-5) v(v-l)+-4 




lfl(v-5)^D-N'-P} 


dD D 


dq (v-5) 


dv 


q(<av-2) dq 


(4.30) 


(4.31) 
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iE 

dt] 


.Dy 


dP 2 
dri ^ 


(4.32) 


dri 


N 


(v-6) 


T- 

( 3 - 0 ) dv 

+ — 

.UJ 

2 dq 


(4.33) 


For the isothermal flow we also assume that the jump condition are heated by isothermal 
flow as given by equation (4. 12-4. 16). 


RESULT AND DISCUSSION 

The transformed jump condition at the shock for both adiabatic and isothermal flow are 


V(l) = 


D(l) = 


f2-\ 




n 


(i-n) 


p(i), 


4v 


0 + 1 0(2 -0)m^ 




Ms 


and 


m 


37 


O 


(1-0)yM^- 11 , [Mi (70-2)0 (2-0)ol 

2 I 


(y-i>i= j ' l2(i-n)(r-i) 

Where V is obtained from the relation 


V = 


(9^ n 
- vM' 

Uj 


2(1+0) 




\^aJ 


(l-£2) 


-1 


(4.34) 

(4.35) 

(4.36) 

(4.37) 


(4.38) 


(4.39) 


For exhibiting the numerical solution, it is convenient to write variable in non 
dimensional form as. 
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iiic iwiuie OI now ana neid vanable tor both the adiabatic and isothermal cases 
are illustrated through figure (4.01-4.05) it is clear fi-om the figure that velocity, density, 
pressure and magnetic field are maximum at the shock fi-ont and decreases rapidly 
towards the centre of explosion, when the flow behind the shock wave is isothermal but 
these variables are minimum at the shock front and increase steadlv so we move townrH® 
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Withoirt magnetic field 
With magnetic field 


Pressure Distribution 
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Without magnetic field 
With magnetic field 


Velocity Distiibution 
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Chapter -V 

Self Similar Solution of Cylindrical Shock 
Wave in Magnetogasdynamics. 

INTRODUCTION 

Self similar solution for shock wave have studied by kopal and Lin |1] using 
similarity method of sedov [2] Elliot [3] wang [4] Helliwell [5] have considered the 
problem of shock wave with thermal radiation. Self similar solution for adiabatic flow 
and is thermal flow was obtained by Rao and Ramana [6] Ojha [7] studied the same 
problem having effects of azimuthal magnetic field. Singh [8] has discussed the self 
similar adiabatic flow of self gravitating gas in ordinary gas. Srivastava [9] is also studied 
self similar solution for magnetic shock wave. G.J. Ball [10] D. Kh ofengein has 
discussed simulation of blast wave propagation over a cylinder & Purohit [11] RF. 
Chishell [12] have obtain an analytic description of converging shock waves D. Bhardwaj 
[13] has been obtain the result formation of shock waves in magnetogas dynamic flows. 
S. Sreekanth, A.S. Nagarajan and S. Venkata Ramana [14] have obtain transient MHD 
free convection flow of on incompressible viscous dissiputive fluid. S. Ahmed and N. 
Ahmed [15] have recently obtain two dimentional MHD oscillatory flow along a 
uniformally moving infinite vertical parlous plate bounded by porous medium. 

Present chapter deals self-similar flow of cylindrical weak shock wave taking 
gravitational forces under non uniform atmosphere taking radiation heat flux into account. 
The radiation pressure and radiation energy have been ignored we have considered the 
problem shock wave reduced by a point explosion in a gas. 

We have considered the problem M. Cheng, K.C. Hung, O.Y. Chong [16] has 
discussed on numerical study of water mitigation affection blast wave Koichi Noguchi, 
Edison Liang [17] has disoissed on three dimensional particle acceleration in 
electromagnetic cylinder and torus. Vishwakarma and Yadav [18] studied the solution for 
self similar flow behind a magnetogasdynamics shock wave in a radiative. 

In this chapter we show that the magnetic field has a significant effect gas is self 
gravitating. The motion of the shock wave is asaimed to satisfy the power law 
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Where A and a are constant and R is the shock radius. The total 
flow increases with time because of the pressure exerted on the gas by 
surface therefore. 


Where p and q are constant and 
distribution law is 


energy. The magnetic field 


Where C and p are constant and the values of q & p are to be determined later, 
The magnetic field is directed tangential to advancing shock front. 

The flow variable just ahead of the shock are 


where 1 + P = a and viscosity is neglected. 

We investigate three types of models, the first having total energy of the explosion 
to be constant {a = 2,5), the second having constant velocity of propogation of shock 

waves (a = 2) and the third having neither constant total energy nor constant velocity 
propagation of shock waves (a = 1.5). An idealized magnetic field is considered for only 
a portion of sphere enclosing the origin i.e. the point of explosion. 


EQUATION OF THE PROBLEM 

The basic differential equation governing. The adiabatic flow in a self gravitating 
gas are given by 
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Where r, t, m, u, p, P, y, 


time, mass contained in a sphere of rj 
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being the shock velocity and ^ is given by the quadratic equation 


The mach number M and Alfen’s mach number Ma are given by ^ and 


Let us seek the solution to the equation in the following form, 


( 5 . 28 ) 

r® t'’ while k, X, a and b are constants and are to be determined from 
the conditions of problem. The total energy E inside the shock wave of radius R is given 


Where 


In terms of the variables q , we can express the total aiergy as 


|(kt5)/a] X-2-b/»(k-5) 


[(2K+5)/a]-I 




+ N - GWOq 


Equation ( 5 . 29 ) yields 
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We choose to be constant at the shock surfaces. This choice fixes the velocity of the 
shock as 


Using equations (5.21 to 5.23) in equation (5.13) to (5.19) and in equation (5.21) we 
assume, without any loss off generality. 


SOLUTION OF EQUATION OF MOTIONS 

The equation (2. 17 to 2. 18) of are used to reduce the equation 


By equation (5. 17 & 5. 18) 
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Uap + u] 


+ O A,t 


putting values u. 


in ^uation we get 


{xa+o: TibJ+u[ko+a'aTi] +n[u 

Using equation (5.32) in above Ration 
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a 


Q' Ti-aOU-SQ' nU+UQ-5U''nQ+OU=0 


5Q'ti ^-uJ + (2-a)uO-5Tiau=0 


50' 


Tif--U|+(2-a)U-5iiU' =0 
il \a J 

Using equation (5.06) 

2!i + i(p+p,)+!k^X=o 

Dt p Sr pr 


(5.34) 


Su Su 1 ap 1 SP^ , h, 5h, h/ ^ 

St SrpSrpSr pSr pr 

u = iu(,,) 


(b) 

(535) 


du. 

dt 

du. 

dt 


= r 


<5t t vt‘ 


t t t^ 


JL 

St "t" 

& 

Su 1 


[bqU'-u] 


[u ap + u] 
Sr t . 

By equation (5.24) 

P = r""-t"-^ P(p) 


(v) 


(ii) 


Sr 


=t 


X-2 


(k + 2)r“.P(n)+r“P(r,)^ 


— = 
Sr 


(k+2)r‘'"'p(p) + r‘‘^lP'(p) 


ap 


t^-2 j.kP! g + 2)p(p)+app’(^ 


(Vi) 


^Pr _ ^k+1 

dr 


fe + 2)P, W+ar|P,'(,,3 


Using equation (5.27) 


h, = r 


_ ^{X- 2 )I 2 

dr 


( k + 2)/2 ^(^- 2)/2 


(n) 


JN,(n)+r ^ -KM^ 


= t 2 


k +2 


k +2 




ari 

r 


X-~2 k 

t 


^ ’N,(q)+N’ (q).aq 


/ k +2 . 


hi = 


r > t = N,{ti) 


J 


putting the value of equation (5.21) & (v) to (ix) in equation (b), we get 
r,, 1 

- jbq U -U1 + 

U . t t ‘ - r‘ t' 


[bn U’ -u]+ -U - [bn u' - u] + ^ Q e-\ [(k + 2)p+an p ] 

t t T \ 


r‘‘''[(k + 2)P^+anP;]-f 


k +2 ^-2 

r 2 t 2 ^ 


+ 


rn 

r 2 ' t “T 

V I 

r'^ t^O.r 


r’' t^O 


X-2 k 

^ t 7 r2 


k + 2 


N, 


[bnu-uj+u[bnu-uj 

+0 [(k +2)P+anP']+ ^ |k + 2)p, +a n Pr^ N, 

substituting the value (5.32) 

k=-a,k=0 b = 4 + q a=-5 & a = 

in above ^nation 
10 


(k + 2) 


N, + N’ 


+N, 


10 
4 + q 


u'n-u + u^-5uu’n + ^^ p- +^— ^ Pr -^^ = o 
a n n a " a 
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(vii) 


(viii) 

(ix) 

+ N'^an 

= 0 

(5.36) 

(5.37) 
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Using equation (5.07) 


By equation (viii) 


Putting these values in equation (c) 


Substituting the values of equation (5.32) 
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K=- a, A,-0, b- 4 + q, a= -5 and a = 


10 

4 + q 


10 ( 2 - 0 .) 


^ 2 U +n4-5U'ti-uJ+N^ U = 0 


5 pn; 


(2 

--U 

Va ) 


■5U'tiN, + (U-i)n, -|UN,=0- 


c K 
5 — ^ T1 

N, 


U 


a 


-5U' n+{u-i)-|u = o 


Now, by using (5.08) 


D u ^ 1. u . 

— hg + hg — — h h« — = 0 

Dt ® “a- % 


^ t, , , , Su hgU 

— he + U + he 1 - - 2 — 

dt 8r ® a- r 


u = ^U(q) 


(d) 

(5.21) 


k.-f2 

h,=r " 


ai 


a 


k+2 

e _j. 2-^ 


X-2 




NeW 


X -2 


+t ‘ N,(ri) 


a 


k+2 


1-2 


X-4 


2,-2 


brj 


2 t ' Ne(ti)+t VNe(Ti) ^ 


Since 


dq _ bq 

a t 


dt 


k+2 


X -^2 


x-4 


X-2 


t 2 N9(q)+t 2 'N0(q)bq 


St 


k+2., , l- 


1-^2 


Ne (q)+bqNe (q) 


k+2 : 2.-2' 


he=r ^ t " N,(t)) 


a 


x^2 


^k+ 2 . ■ 

r 2 Ng (q) ^ (^l) 

r jL 


(xi) 
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Substituting the values of equation (5.21 ), (5.26) (ii), (xi) (xii) in equation (d), we get 


We get 


SU'riNe +(U-1 )n, 


Now, using equation (5.09) 


By equation (5.24) 
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p = P (p) 


^ _ j.k+2 

dt ^ ■ 


a ^ ^ 


_k+2 


t^-2 p' ^ +(x_2)pt 


X-3 


5t , 


[(^-2)P + P'bp] 


Again by equation (5.24) 
P = r‘“ P (t,) 


» = ,« 
a 


P' ^ +P(k + 2)r ’'^' 
dr ^ ^ 


a-2 




,k+l 


dP 


^X-2 j.k+1 


dr 

By equation (iv) 


P aq +(k + 2)p] 


ap t^r" 


dr r 


_kQ +0 aq] 


By equation (v) 

Putting values of equation (iv), (v) (xiii), (xiv) in equation (c) 
t^'" [(X-2)P+P’bq]-y - 




jfi bq+A,oJ 


^k+2 jp 


rt"Q(q) 

+ [p' aq + (k + 2)p] 


j.k+2 jX-2p ^l j.k 


y L£_ [k O + O' 

aq] 


(xiii) 


(xiv) 


(iv) 


(V) 


[(X - 2) p + P' bq] - y [□' bq + Xo] ^ + ^k + 2) P + P' aqj -y [k O + Q’ a q] — = 0 

Q O 


Now using equation (5.32) 


lU 



We get, 


Now using equation (5.10) 




Putting these values is equation (f) 
t^ |(k + 3)W+ W' aTiJ= 47E. r 
(k+3)w + W’aTi = 47ta 

(k + 2)w + W’ aT|-^43iQ = 0 


The boundary condition (5. 12) to (5.18) are transformed into the following form 
Boundary condition (5.12) is 



Boundary Condition (5. 13) is 
Pi =Po^ 

r'^ t^ Q (ti) = AR““ ^ 
r = R, K= -a, X, = 0, r\ = tig 
R'“ t° O (t1o)= AR'“ 5 







WiMi 


Boundary condition (5.14) 


p(n) ='P AR 


t^-'p(Ti)='P AR 


'PAR 


Boundary condition (5. 1 5) 
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Boundary Condition (5.17) 


Boundary condition (5.18) 
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w(t1o) ^ 4 Tt 

A (3 - a) 


w (n.) 


4n 


( 3 -a) 


RESULT AND DISCUSSION 

Similarity solution of the problem of propagation of a cylindrical 
magnetogasdynamic shock wave have been obtained. For numerical calculations the flow 
and field variables have been taken in following non dimensional form. 


u 


r| 


U, 



•N 

_p_ ^ 


Pi 

UoJ 


r 

_L = 

ti 

P, 




_p_ ^ 

r] 

P. 


^6 

n 


UoJ 

ik - 

f \ 

n 

h 

Zj 

UoJ 


m 



-a /a 




inoj 


u(ri) 

u(rio) 

Q(ti) 

“ p(ti) 

Pr (^oj 
Pr(n) 

Ne(Ti) 

Ne(no) 

N.(ri) 

W(p) 

wk) 


The differential equation 3.04 to 4.0 are numerically integrated by well known 
software matlab and computer analysis made for flow variables through graphs for certain 
choice of parameter y =1.4 = 20, =10 & a = 0.25. It is found that variation in 

velocity, pressure density are minimum of shock fi'ont but increase towards the line of 
i ocplosions. But variable magnetic field of mass variable is maximum in shock fhwrt but 
decrease toward line of explosion. 
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Chapter -VI 

Analytical solution of spherical shock wave in a 
rotating gas with axial con^onent of magnetic field 

INTRODUCTION 

The propagation of hydromagnetic shock through a rotating gas being close to the 
actual situation of thunder. Kumar and Prakash [1,2,3] recently using C.C.W. [4,5,6] 
method have investigated the propagation of weak and strong diverging cylindrical 
hydromagnetic gas taking desntiy distribution variable but axial component of magnetic 
field of constant strength. 

In this chapter, using the C.C.W. method we have considered the problems of the 
propagation of a converging and diverging spherical shock wave through a rotating gas 
under the influence of a magnetic field of constant axial, simultaneously for both weak 
and strong waves we have assumed initial distribution is constant. Sharma V.D. [7] & et 
al. have studied the problem of strong shocks in an ideal gas [7], [8], Pandey [9] use an 
approach based on interaction of an characteris tic shock with a weak discontinuity in a 
non ideal gas. T. Raja [10] also find similarity solution, for waves. 

The analytical expression for shock velocity and shock strength have been 
obtain^ for weak shock under the two conditions namely when the magnetic field is 
weak and when the magnetic field is strong. For strong shock also we have considered 
two cases, i.e. when the magnetic field is strong and when / /?, is approximately equal 
to f + 1 / - 1 which is purely a non magnetic case. 

BASIC EQUATIONS BOUNDARY CONDITION AND AN ANALYTIC 
EXPRESSION FOR SHOCK VELOCITY 

The equation governing the spherical symmetrical flow of rotating gas under the 
influence of axial component of magnetic fields are written as 
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5u 5u h 0p 

_ + U + — -E- 

dt dr p dr 


(d dp^ f . 

W = o 


dp dp ( du 2u 


v‘ u. 5 , 

~ ~ ~ = 0 
r 2p cr 


+ u 1 - p 

dt dr ^ 


dr 


+ 


( 6 . 01 ) 


( 6 . 02 ) 


(6.03) 


dH, dH. 


+ 


dt dr 


+ H, 


( du u'] 

dr r 


(6.04) 


where u, p and p are the velocity, pressure and density behind the shock front 
respectively while, is axial components of field H. 

The magnetogasdynamics shock conditions can be written in terms of a single 
parameter. 

N = p, /po as 


p, = PoN, H, =HoN, U, = 


N- 1 

N 


U 




2N 


(t+i)-(y-i)n 


ao + ^+{(2-y)N+y} 


p =p ^Po(N-0 
' (y + l)-(y-l)N 


2 . y-1 


ao + 


hi {N-if 


(6.05) 


where 0 and 1 denote respectively the states immediately ahead and behind the shock 
front , U the shock velocity, a the sound speed (ypo /pg) and the Alfven speed 


IS 


where 
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In the equilibrium state we have 


JL _ X! 

Po dr r 


= 0 


Integrating the proceeding equation we get. 


P = K] log — and a^, = K, 

r 


where K, = p' 


log 


K, 


P' . 


and 


p' V*- log A is constant of integration. 

WEAK SHOCKS 

For a very weak shock we take the parameter N as 


N= 1+ s, 


£ <<1 


Now consider the two cases of weak and strong magnetic fields. 
Case -I 

When the magnetic field is weak i.e. 
bo «al (i.e.pH^ «TPo) 
and therefore, the boundary conditions (5) reduce to 
P, = Po (i + e) 

Ui = bo • 


132 


Chapter- VI 


Hzj =Hzo (l + s) 

P1 = Po(1 + t) 


f 


u 


y + 1 


1 + 

V 4 


(6.06) 


Case 2 


When magnetic field is strong » ag under these circumstances the 
boundary condition (6.05) reduce to 


Pi 


■ Po 0 + s) 


U, = bg 


Hz, =Hzo (l + e) 


Pi 


U 


Po 0 + s) 


3 

1+-8 , 

4 ) 


(6.07) 


STRONG SHOCK 
Case -I 


When the magnetic field is strong b^ » a^ (i.e. p »YPg) under this 
condition the boundary condition (6.05) becomes. 

Pi =Po N 


H, = Hg N 


U, 


N-1 

N 


U 


Pi XU' 


Po 


+ L 
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where 


X 


Y (Y+l)(N-iy 

2N{(2-y)N + y} 

(y+i)n-(y-i) 
(y + 1)-(y-1)N 


and 


and 


(6.08) 


Case -II 


For strong shock when (i.e. when N -»[(y + i)/(y - 1)]) is small the field 

may be r^arded as independent of magnetic field. 

For converging shock the characteristics form of system of equation (6.02) is 


dp + p (Hj, dHj,) + pcdu + 


pc^u ^ pcv^ 


u - c u - c 


dr 


0 


where as the characteristic form of diverging shock is 


dp + p (h^ dH^)+ pcdu + 


pc^u pev^ 
u+c u+c 


dr „ , 
— = Od 
r 


(6.09) 


( 6 . 10 ) 


where c^ = +b^ = ^ + iiS 

P P 


and 


WEAK SHOCK WITH WEAK MAGNETIC FIELD 

Substituting shock condition (6.06) into (6.09) and neglecting the second and 
higher order term of s since e < < we get. 




d + 


dpo dao dr 2dpo 


Po 


YPo 


= 0 


( 6 . 11 ) 
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Now substituting 


i!?£, a«d 

Po ^0 YPo 


in equation (6. 1 1) we get 


ds 

& 


[y-4] , 1 


rP? Pf 


+ — log 


dr 


( 6 . 12 ) 


where P^ = 


pH 


zO 


yk, 


pf = o 


Integrating of equation ( 6 . 12) yields 


8 = K 3 rf^-« exp 


log 


2P( ^ rj 


(6.13a) 


where log to constant of integration in case of diverging wave shock with weak magnetic 
field, we have the relation. 


s = K 3 r 


- 1/2 


A 

log- 

r 


V 'll 


exp j- I Pi (log)—)'’ 

8 r 


(6.13b) 


where log K 3 is constant of integration and 


Qi = 


4 4 




Hence from equation (6.13a) and (6.13b) we have 


2 L = 1 + I±> K, rW-dWiR- exp. 


1 (, A' 


(6.14) 


I±iK3 r- 

ao 4 ^ 


r aV"‘ [ 3 r, aV 

^log- expj- -pf log- 
T J 0 V t/ 


(6.15) 


for converging and diverging shock respectively. 

WEAK SHOCK WITH STRONG MAGNETIC FIELD 

Substituting the shock conditions (6.07) in (6.09) we get 
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ds + 


dpo _ 

, Po ypo 


2 f 


dr db„ 


+ ' 


V 


2dpo 


= 0 


0 ) yPo 


(6.16) 


Substituting the respective quantities in equation (6. 16) we get 


8 


Y(l + Pf- 2 ) 

y 


^ A 
log - 




r 


(6.17) 


integrating equation (6.17) we get 


e = K. 


log — 

V r j 


-S2 


(6.18a) 


where log is constant of integration and 
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lr(i+p;)-2l 


when we consider diverging shock with strong magnetic field we get 


8 = log 






■1 exp 


2pf (l-pf)flog 


(6.18b) 


where log is constant of integration equation (6. 1 8a) and (6. 1 8b) disclose that 


u 

b„ 


1 + -K. 


log 




(6.19) 


u 3 ^ Yi 

— =l+-K^ log — 

ao 4 V rj 


-lexp ]-2p," (l-pj^Hog 


A 

r 


( 6 . 20 ) 


and diverging shocks respectively. 


Chapter-Vl 


STRONG SHOCK 

Substituting shock condition ( 6 . 08 ) into ( 6 . 09 ) we get 



dU' 

X N-l 

fx^ 

1/2 i 

+ U" 

'x dpo 

2x dao ^ x(N- 1) dri 


_y 2 



.Y Po 

y ao (n-i)-V^ rj 


+ 


dpo yk, 


Po 


Po 


+ 




| (N-1)+/'nx ] dr , ^ N-^N^ 

■1)-/'Nx 


= 0 


Substituting the respective values we get. 


(6.21) 



dr r 
where B = B’/M 
C = C’/M 

p,_ x(n-i) 
(n-i)./” xn 


( 6 . 22 ) 


M 


(N-l) 




u. 


and 


^,_Lk, ^ yk^N^pf f (N-l)+-/'NX ] n/nX 

p' p' |(n-i)-V~nxJ (n-i)-/'nx 

Integration of (6.22) yields 

U^-K^r-P+I (6.23) 

where K j is constant of integration. In the case of diverging strong shock with strong 
magnetic field we have 
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U2 = + 


B (6-24) 

Where K 5 is constant of integration, finally for converging and diverging shock we have 


1 




K^log- 

r 


B 


(6,25) 


and 


U' 


1 


log - 


K, r'P + 


-1 

B. 


(6.26) 


respectively 

In the last, the expressions for the velocity, the density, and the pressure of the gas just 
behind the shock surface are 


STRONG SHOCK WITH WEAK MAGNETIC FIELD 


U = K,K 3 log 

" V r ; 


^ AY ^ (r-4)/TBf J 1 , TaY 

’ \t J 


2 P? 


(6.27a) 


U = K 2 K 3 


1/2 


log — 
V rj 


exp |-^P? log|^yj 


(6.27b) 


P =P' 


I + K 3 exp 


Vflogf 


Vi 


2 p^ 


(6.28a) 


P =P’ 


I + K3 r'*'^ 


r, 

r 

3 „2 

f, A1 


log — 

exp i 


log — 



1 

8 

1 r J 

J 


(6.28b) 


p=k, 


log— + ykj 
r 


V 


fl 

log — 

V ry 




exp 


2Pf 


log 


(6.29a) 
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(6.29b) 


(6.30a) 


(6.31a) 


(6.31b) 


(6.32a) 


(6.32b) 


(6.33b) 


p = p' N for converging diverging cases 
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Chapter-VII 

Self Similar Model of Radiative Shock Waves in 
Magnetohydrodynamics with Magnetic Effect 


INTRODUCTION 

In this chapter we consider propagation of shock wave under the effect of 
radiation taking magnetic effect in to account of an instantanous release of energy in a 
non-uniform equilibrium conditions. The disturbance are headed by a shock of variable 
strength. This model is of considerable physical interest in sonic booms. 

The medium ahead of the shock is assumed to be inhomogeneous and at rest the 
shock position in this problem is given by 

The theory of shock waves and related flows has attracted a renewed interest in 
connection with the phenomena associated with laser production of plasma, nuclear 
detonation and astrophysical situations. Ray [1] obtained an analytic solution in the case 
of a central explosion in a gravitating mass of equilibrium in which the disturbance was 
headed by a shock surface of variable strength. Ojha [2] considered the propagation of 
explosion waves in a steller model in the shock strength does not remain constant in 
general and mach number of the shock is a function of time. Ojha and onkar [3] studied 
the propagation of shock in inhomogeneous self gravitating gaseous mass are headed. 
Tayar [4], Carrusetal [5], Klynch [6], Sedov [7] have accumulated and extensive to 
literature on self similar model of phenomena for the propagation of shock wave in gas 
dynamics and in the formation of stars and super nova explosions. Summers and whitham 
[8] have discussed qualitative analysis of self similar solution to the problem of unsteady 
spherical symmetric motion of self gravitating gas Sakurai [9] has considered the problem 
of shock wave arriving at the edge of a gas in a medium in which density varying as 
power law. The flare generating shock waves has been studied by Hundhanses [10], Deb 
Ray and Bhownick [11] have obtained the self similar solution for central explosion in 
stars with radiation flux when the shock is isothermal & transparent. 
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where p and A are constant & p < i 

We assume that the density disturbance is given by 

Po = (7 02) 

where b& P are constants. 

The magnetic field disturbance the magnetic field distribution is taken cc & 
Rosenate & Frankenthal [12] as 


directed tangentially to the advancing shock fi-ont where c & 5 are constants. 


EQUATIONS OF MOTION 


Where p, u, p, hg, h^, y,m are the density, velocity, pressure, azumethal 


magnetic field, transverse magnetic field and ratio of specific heats of the gas, mass 
respectively. 


We introduce the similarity variable t] 


and take the solution of 


fundamental equation in the form 
u=R V(ri) 
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he = R Hq (ti) 

h,=V^RH,{Ti) 


(7.12) 

(7.13) 


dR 

where F, D,P,Hg,H^ are function of p only and R = — the shock velocity. 
The jump condition for a shock wave (cf summer & Whitworth [13]) are 


u, 


P, = 


2R 

(y+l) 

2PoR^ 
(y + l) 


_ (y + i)po 

A = 

(/-l) 

h 

(r-i) 


(7.14) 

(7.15) 

(7.16) 

(7.17) 

(7.18) 


where suffix 1 denotes the volumes of flow variables immediately behind the shock front. 


SOLUTION OF EQUATION OF MOTION 

By Equation (7.01) 

dp dp pdu - 

dt dr dr 

By condition (7.09) 

u = RV(ti) 


^ -R V' 


(7.01) 




dr 


dr 


' R(t) 

Sp 1 
___ 
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By condition (7. 10) 


By condition (7.10) 


By condition (1.02) 


substituting these values in equation (7.01) 


R + VR 
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V 

-Of. + — 

i? R 


D' + 


Dp 

n 


DVt] 

+ =0 

rjR 


■Dr} + V 


D + 


n 


+Z)F' =0 


{-r} + v)D 


''±.v 

n 


D = 0 


(7.19) 


By equation (7.05) 


dp 


dP yP f dp 


+ u 

5t 5r p 
By Condition (7. 1 1) 
P=Po P(ti) 

DP 


+ u 


dp 

dr 


0 


^ ^ Al 
dr R 


D +■ 


n 


substituting these values in equation (7.05), we get 


R 


-Pr}+2P 




kvp,ie 


R 


-p 


n . 


rp,Fep p^DijR VRp, 

• i -' 


AO 


R 


R 


D + 


DP 




=0 


{v--n)p' +2P^y-i + 




■pv ^ 2(p-l)' 

p-^r 


D 


7} 


(7.20) 


Now, taking equation (7.06) 
5u 


du 1 SP (he + h ) d / , \ 

a dr p dr pr a ^ ® ^ 


By condition (7.09) 

u = RV(ri) 

da _ RV 
dr ~ R 

4^ 


dt 


P' iP + 2 R RP 

dt 




-r]R 
V ^ y 


R^^IRRP 



By condition (7.01) 


again, by condition (7. 1 1) 


By condition (7.12) 

K =4^ RHg (tj) 
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again, by condition (7.13) 


By condition (7. 11) 
P = p, R^P(il) 


substituting these values in this equation 


PoD.qR 


Using equation (7.07) 
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H,^R + HgR 
dt 


PQ 




n 


dH, 


R 

ie 


R.R+ H, 


{p-^W 


juR 




P 


again, h^ = ^ R (t|) 


dhc 


dr 


R 


H 


5ti 


dr 


e “ -vPo + Hq 


SPfl 


dr 


R 


H, 


Po . Hep 


+ 


R 2tiR 


Po 


dhg _ R 
dr R 




HeP 


2 // 


By condition (7.09) 
u = R V(tj) 


dr dr 


dr 


R 


substituting these values in eq. (7.07) we get 






H.+ 


HgP 


+ 


^ RH^ 
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^ + H H,=0 

^ /j 2 ri 

By equation (7.08) 

dh. dh da 

— L + u — ^ + h_ — =0 

8t dr dr 

By condition (7.13) 

h =^R H, (n) 


(7.22) 


^=Jpo H^^R + H,R 
dr ^ ' di 


^ " R ' q R 


dt R ‘ /i ^ 


lin h, = ^ R H, (t]) 


dr cr dr 


= R 


dh, _ R4 p^ ff' ,iM. 

dr R ^ 2ri 


By equation (7.09) 
u = R V (q) 


^ = kv' 2 IL 
dr &• 


du _ kv 

dr R 


Substituting these values in equation (7.08) 
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Equation (7. 19-7.22) can be reduce in new form 


P + yPV 


(V-q) y( 2 + |3)V pV 


piV-qf-yP+^He+HJ J (7.28) 

Now shock condition 7.14-7.18 are transformed into following form by condition (7.09) 
and (7.14) 


151 


Chapter- Vn 


v(i) = 


(y + l) 

By condition (7.11) & (7.15) 

2 Po 


(7.29) 


P, = 


(y + l) 


Po 


R^P(tl): 


2PoR^ 

(y + l) 


at T\ = 1 


P(l) = 


(y + l) 


(7.30) 


When the magnetic field is weak J.B. Singh and S.K. Pandey [13]. 

^ (y + 0 
Po (t-i) 

Which is purely non magnetic By condition (7.10), (7.16) and (7.31) 
(y + l) 

at p = 1 

D(l)=l 

By condition (7. 12) and (7. 17), (7.3 1) 

(y-1) ” 

at p = 1 

He(l)=I 

By condition (7. 13), (7. 18) & (7.3 1) 

(y+0 

(y-i) 

p = i 

H.(l)=l 


(7.31) 


(7.32) 


(7.33) 


h = h. 


zO 


(7.34) 
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And appropriate shock condition are 
2 


V(l)^ 

P(l): 


(y+i) 

2 


(y+l) 
D(l)=l 

He (l)=l 


(7.29) 

(7.31) 

(7.32) 

(7.33) 

(7.34) 


RESULT AND DISCUSSION 

The numerical result for certain choice of parameters are in tabular form. The 
nature of field variables is in illustrated through tables for following sets of parameter. 
7=4/3, 7/5, b = 3, p, = 0.7 p=0.4. It is obvious from table that velocity, density, and 
magnetic fields increase as we move towards the point of explosion while the pressure 
decrease as we approach towards the centre. The variation in flow variable is negligible 
for different value of y . 
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0.995 


0.994 


0.992 


Table-7.2 























mfl 

■I 


0.994 

0.993 


0.992 


0.991 


0.747810 


0.752018 


0.755959 


0.759653 

0.763123 


0.766386 


0.769459 


0.945030 


0.957211 

0^967938 


0.977202 

0.984998 


0.991333 


0.996217 


Table-7.4 


0.972216 


0.929031 

'0.985778 


0.942585 


0.999573 


0.9568/55 


0.914534 





0.927960 


0.955445 


0.982419 


0.908850 

0.934714 


0.959989 


0.984661 


0.908721 


0.932167 


n 

V 

D 

i 

1 P 


1.000 

0.733333 

0.900000 

0.900000 

0.900000 

0.999 

0.742203 

0.940162 

0.933239 

0.952128 

0.998 

0.749851 

0.975842 

0.966437 

0.901436 

0.997 

0.756477 

0.906988 

0.900357 

0.947808 

0.996 

0.762249 

0.933637 

0.935645 

0.991201 

0.995 

0.767304 

0.955903 

0.972824 

0.931637 

0.994 

0.771756 

0.973957 

0.912297 

0.969188 

0.993 

0.775699 

0.988014 

0.954354 

0.903970 

0.992 

0.779212 ^ 

0.998318 

0.999186 

0.936128 

0.991 

0.782362 

0.905131 , 

0.946894 

0.965828 
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Chapter -VIII 

Similarity Solution of Cylindrical shock wave with 
radiation energy & material pressure in 
Magnetohydrodynamics 


. r 


' 


INTRODUCTION 

The theory of shock wave and related flows are considerable physical interest 
Shock waves conceivably driven by solar flares are observed to propagate into 
interplanetary medium. Rogers [1], Deb Ray [2] have obtained an exact analytic solution 
for the shock wave problem with an atmosphere of varying density. They have considered 
the problem taking effect of radiation energy and material pressure shock propagation at 
very high temperature in which the radiation effects might play a very important role 
through the coupling of radiation and magnetogasdynamics fields on account of the high 
temperature, gases are ionised over the entire region of interest in the shock and the 
medium behaves like a medium of very high electrical conductivity. The explosion along 
a line in a gas cloud has been discussed 

The propagation of shock waves has been studied by Greifigner and Cole [3] and 
Green span [4], Christer [5] and Ranga Rao & Ramana [6] without taking into account the 
radiation effects. Elliot [7], wang [8] and Helliwell [9], have considered the effect of 
thermal radiation in their studies of gas dynamic using similarly method of sedov [10] 
Singh [11] has studied the problem in ordinary gas dynamics. Theoretical and 
Experimental studi^l of radiative shock have been treated by C. Michaut & et al. [12]. 

In this chapter we consider the problem of cylindrical shock wave in 
magnetogasdynamics when the atmosphere is ncm-uniftMin and conducting taking counter 
pressure and radiation flux into account. The radiation pressure and radiation energy have 
been considered. The gas in the undisturbance field is assumed to be at rest and it is grey 
and opaque. 

We suppose that the magnetic field Hq and density p, distributions ahead of 
shock vary as an inverse power of radial distance from centre of symmetry i.e. 


/■ 


iliiiilf; 


; ‘i ' ' ' / 
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«o =- 

6, j. 


and 


B 


(-2<w<2) 


where A, B and w are constant 


SELF SIMILAR FORMULATION 

The cylindrical polar coordinate, where r is the radial distance from axis of 
symmatry are used here. The equation of conservation of mass, momentum, energy and 
magnetic flux in the infinite conduction region behind the wave are. 


dp P ^ t \ n 

— + — — (ru) =0 

dt T dr ^ 

dt dr r a- ^ 


^ + He — = 0 

dt ® dr 


+ H .^=0 

dt “ ar 

d { \ ( \ d ( \^ 


( 8 . 01 ) 




du d ( \ 


KP) 

pHe d 

r dr 


+ - I: (r^) = o 

pr Sr 


( 8 . 02 ) 

(8.03) 

(8.04) 




where 

d d 


dt 


d 

— + u — ■ 
a dr 


(8.05) 


and p is the density p the pressure, u the radial velocity, h the azimuthal magnetic field; q 
the heat flux, t the time and e is the material energy radiation energy pr radiation 
pressure. The magnetic permeability p is taken to be unity. For an ideal gas we have 

PP 


(y-i)’ 


FpT 


(8.06) 
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Where y is the adiabatic gas index, T the temperature and y the gas constant. Assuming 
local thermodynamic equilibrium and taking Rosseland is diffusion approximation we 


is Stefan Boltzmann constant; C the velocity of light; and v, the mean free path 


where 


of radiation a function of density and temperature following wang [13], we take 


where , a and p are constant 


Given P 


By condition (8.07) 


Substitute the value T 


substitute the value of ‘v’ in above equation using 


Where 


Following Singh & Vishwaka [14] the jump conditions at an isothermal shock front are 
taken to be 
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where 


+ 2/ M a 


SIMILARITY SOLUTION 


By the standard dimensional analysis of Sedov [10], the non dimensional variables 
T| is defined by 


where S 
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and the dimensionless constant m is defined such that f] assumes the 


value one on the shock front. This choice enables us to write ti 


where R is shock radius. 

Then the field variables of the flow variables of the flow pattern in terms of 
dimensionless function of -p are 


SOLUTIONS OF EQUATION OF MOTION 
By equation (8.01) 


By equation (8.19) 


By equation (8.22) 
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By equation (8.19) 


putting these values in equation (a), we get 


IG-G' 5ri+UGri-4GU+UG+U'Grj + 0(7 = 0 
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G' T]iU-S)-2GU + 2G + U'Grj 


dG 


G 


2(U-l)- 


11 


dr| 


dp p(U-6) 

Now by equation (8.02) 

^ dt 8v T dr ^ 


(8.25) 


da da dP 5H 

p + pu_ + -— + pHg + 

ot or 0r dr 

By condition (8.21) 

A 


e , 


H 


He (n) 


(b) 


dr 

dr 

dH. 


A , dll 

a 2 

r dr r 


M ’ ^ ^ f4 

n g — -Y n g 


A 

.2 


dr r 
By equation (8.20) 
^ A" 


[H.rt-H,] 


P(p) 


dr 


r dr r 


Z_ H 

2 

r ' r 


Sr 

& 




[p’77-2P 


du 1 




By condition (8. 19) 


-u ^ u_ 

t t 


du r 


dt t 


^ [^t/' Sr,-u] 


substituting these values in equation (b) 


A't' „ r 


G -L [_u’ 6 ti-u] + G. - U 


G[-U-U'5q + U'+UU'qJ+ [p'ti-2p]+ ^iHe H 

dr] (t/-<5) [g ^ \ dr] r] iG 

By equation (8.03) 


dt ® a- 


SHg 0Hg 5u „ 

+ u — - + Hq — = 0 

dt dr \dr 

By condition (8.21) 

Hq = — Hg (q) 


_r, w. an 1 ... 
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By condition (8. 19) 


\j + U'ri 


Putting these values in equation (c) 


H^A U + U'fj 


By equation (8.04) we get 
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By equation (8.06) 


By equation (8.24) 


By equation (8.22) 


By equation (8.22) 




-GP, +GP; tj+G'P 


By equation (8.23) 


substituting these values in equation (d) we get 


[2(f+f, )g-( 


£/. [-g(/>+/>, )+G (;>’+/>;); + c (/) 


Now by equation ( 8 . 05 ) 
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du d 

P + — — 

dt dr 


(/>+/>, ) 


+ 


d 


»^K)=o 


dll dP ^ dP^ HH. dH 


P __ j ^ 

dt dr dr 


+ 


e , pHl 


dr 


+ 


^ Su ^ ap 5Pr „ aHe fiH' 

St dr ar ar r 


By condition (8.19) 


u = '-U (v) 


t 


du r 

dt p 


\-U' 5t]-u\ 


By condition (8.20) 
A 


P = ^P(n) 

[p't]-2p\ 

By condition (8.24) 


^ =i 

dr r' 




dP 


=il 

dr 


^'rV- 2 Pr] 


By condition (8.21) 


He = T He W 


dH, 


■dr 


A 

^2 


\h’s t] - Fg] 


substituting these values in equation (e), we get 




G 


r 

77’ 






(e) 


A 2 ' .ry^ 

+ 11 - 3 - He- =0 

r r 






ill 
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Putting of values p =1 and H' 


where k 


a non dimensional radiation parameter. 

The transformed shock conditions are, using equation (8. 1 1-8. 16) & (8. 19-8.24) & 
(i), (ii) By condition (8. 14) & (8. 19) 


By condition (i), (8.13) & (8.21) 


h.(0=n 

By condition (8. 1 1) & (8,22) 
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By condition (8.12) & (8.20) 


By condition (8.15) & (8.24) 


N. M fuH 


By condition (8.16) & (8.23) 


m 
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where 


as the initial values for our numerical calculation 


RESULTS AND DISCUSSION 

In order to exhibit the numerical solution it is convenient to write the field 
variables in the non dimensional form as 
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U/U2 

P/P2 

1.200000 

1.200000 

1.189713 

1.211770 

1.179343 

1.237695 

1.168883 

1.258038 

1.158326 

1.283072 

1.147663 

1.303076 

1.136886 

1.328336 

1.125985 

1.349134 

1.114949 

1.355754 

1.103765 

1.378466 

1.192420 

' 1.397525 

1.180898 

1.423162 


1.16918 

1.1572S 

1.145091 



1.491199 


Table-8.01 

K = 10,6=2 

P/P2 

1.200000 


1.233358 
1.252275 ’ 
\ 211232 
1.298743 
1.317355 


1.343649 

1.388237 

1.391761 

1.424887 


1.454297 


2.512686 


2.538742 


2.557137 


2.578501 I 


H/H2 

1.200000 


1.259207 

1,288322 

1.317315 


1.346149 

1.374785 

1.413179 


1.431281 


1.459036 


1.496383 


1.523254 


1.549573 


1.575256 


1.590211 


1.614334 




2.205242 

3.225350 

4.252065 


7 


7.312421 

8.339657 


10.350662 


11.377211 


13.491037 


15.513795 


17.537041 


18.55182 


20.570432 


22.612750 


1.201495 

1.104061 

1.207675 


1.212310 


1.217943 

1.224552 


1.232112 


1.240601 

1.249999 


1.271435 


1.209902 
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Tabie-8-02 

K = 10,5=5 


n 

U/U2 

P/P2 

P/P2 

H/H2 

q/qj 

T/Ti 1 

1.00 

1.200000 

1.200000 

1.200000 

1.200000 

1.2CX)CX>0 

l.”20CKX)0 

0.99 

1.180721 

1.227090 

1.308584 

i. 402352 

i.557591 

1.551374 

0.98 

1.171376 

1.252716 

1 356953 

1.424803 

2.575334 

1.593582 

0.97 

1.161957 

1.277757 

1.456264 

1.457320 

4.597607 

1.616606 

0.96 

1.152454 

1.303175 

1.487805 

1.479862 

5.619113 

1.65M27 

0.95 

1.132855 

1.350017 

1.493009 

1.502382 

6.664906 

1.705027 

0.94 

1.093147 

1.459420 

1.534625 

1.534821 

7.690399 

1.720391 

0.93 

1.083317 

1.492608 

1.550913 

1.637111 

9.701375 

1.766503 

0.92 

1.073348 

1.530890 

2.597283 

1.699173 

10.733997 

1.793348 

0.91 

1.053221 

2.585648 

2.614659 

1.710912 

12.754799 

1.820913 

0.90 

1.042917 

2.628327 

2.655302 

1.752220 

13.790675 

1.859184 

0.89 

1.032411 

2.666048 

2.701621 

1.782970 

15.808854 

1.878150 

0.88 

1.021676 

2.693373 

3.736155 

1.793013 

17.826841 

1.897801 

0.87 

1.010680 

3.718656 

3.751529 

1.822180 

19.852347 

1.908128 

0.86 

1.009388 

3.737572 

3.790382 

1.880276 

21.893165 

. 

1.929125 
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